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ABSTRACT 

From the list of 2321 transiting planet candidates announced by the Kepler Mission, we select 
seven targets with favorable properties for the capacity to dynamically maintain an exomoon and 
present a detectable signal. These seven candidates were identified through our automatic target 
selection (TSA) algorithm and target selection prioritization (TSP) filtering, whereby we excluded 
systems exhibiting significant time-correlated noise and focussed on those with a single transiting 
planet candidate of radius less than 6 i?® . We find no compelling evidence for an exomoon around 
any of the seven KOIs but constrain the satellite-to-planet mass ratios for each. For four of the seven 
KOIs, we estimate a 95% upper quantile of Ms/Mp < 0.04, which given the radii of the candidates, 
likely probes down to sub-Earth masses. We also derive precise transit times and durations for each 
candidate and find no evidence for dynamical variations in any of the KOIs. With just a few systems 
analyzed thus far in the on-going HEK project, projections on rjf^ would be premature, but a high 
frequency of large moons around Super-Earths /Mini-Neptunes would appear to be incommensurable 
with our results so far. 

Subject headings: techniques: photometric — planetary systems — stars: individual (KIC-11623629, 
KIC-11622600, KIC-10810838, KIC-5966322, KIC-9965439, KIC-7761545, KIC- 
11297236; KOI-365, KOL1876, KOL174, KOL303, KOI-722, KOI-1472, KOI-1857) 



1. INTRODUCTION 

The "Hunt for Exomoons with Kepler" (HEK) project 
is the first systematic survey fo r moons around plan - 
ets outside of our solar system (jKipping et al.l [2012bl ). 
Whilst many planets around our Sun host one or more 
satellites, there is no empirical evidence for moons 
around the hundreds of extrasolar planets detected in 
recent years. At best, one can interp ret the possible 
detect ion of a circumplanetary disc by iMamaiek et al.l 
as a putative moon- fo rming region. The Kepler 
Mission (jBorucki et al.l l2009f ) is the most suitable in- 
strument available for detecting exomoons thanks to the 
large number of target stars, long temporal baselines, 
nearly continuous monitoring and very precise photom- 
etry. By monitorin g the timing of exoplanet transits, 
IKipping et al.l (|2009D have estimated that Kepler should 
be sensitive to ^^[Mq] mass exomoons. In addition, Ke- 
pler is designed to detect 0[i?©] radius transitin g bod- 
ies and moons may be found in a similar way (Kipping- 
l2011a[) . It is therefore argued that Earth-mass/radius 

^ Harvard-Smithsonian Center for Astrophysics, Cambridge, 
MA 02138, USA; emaih dkippingQcfa.harvard.edu 
2 NASA Carl Sagan Fellow 

^ Dept. of Astrophysical Sciences, Princeton University, 
Princeton, NJ 05844, USA 

* Niels Bohr Institute, University of Copenhagen, DK-2100, 
Copenhagen, Denmark 

° Centre for Star and Planet Formation, Natural History Mu- 
seum of Denmark, University of Copenhagen, DK-1350, Copen- 
hagen, Denmark 

° Citizen Science 
Alfred P. Sloan Fellow 

* Packard Fellow 

® Dept. of Space Studies, Southwest Research Institute, 1050 
Walnut St., Suite 300, Boulder, CO 80302, USA 
t Based on archival data of the Kepler telescope. 



moons should be detectable. Although there are no 
moons this large or massive in our solar system, HEK 
seeks to answer whether this is true for all exoplanetary 
systems or not. 

A detailed description of th e goals and methods o f the 
HEK project are discussed in IKipping et al.l (|2012bl ). To 
date, the analysis of only one system fo r exomoons has 
been published by the HEK project in iNesvornv et al.l 
(|2012) . In this case, the target planetary candidate, KOI- 
872.01, was identified as being a target of opportunity 
(TSO) due to the presence of very large transit timing 
variations (TTV), enabling us to detect a second non- 
transiting planet in the system (and confirm the plan- 
etary nature of KOI-872.01). The non-moon origin of 
these TTVs demonstrated the importance of the careful 
interpretation of dynamical effects. 

In this work, we will provide an analysis of seven plane- 
tary candidates (KOIs) identified through our automatic 
target selection (TSA) method. Each candidate therefore 
satisfies the criteria of having the capability to host an 
Earth-mass moon plus sufficient signal-to-noise to make 
such a detection feasible. Consequently, null-detections 
have much greater significance for understanding the fre- 
quency of large moons around viable planet hosts, "qq^ . 

2. TARGET SELECTION 

2.1. Automatic Target Selection (TSA) 

2.1.1. Overview 

The HEK project treats the Kepler Objects of Inter- 
est (KOIs) as a list of potential moon-hosting targets in 
much the same way that Kepler itself treats the Kepler 
Input Catalogue (KIC) stars as a list of potential planet- 
hosting targets. At the time of writing, 2321 KOIs have 
been reported bv iBatalha et al., (,2012.) (B12). Search- 
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ing individual systems for signs of an exomoon is a time 
expensive task in t erms of computationa l demands and 
human manpower (jNesvornv et al.ll2012f) . As a result, 
the HEK project performs a target selection (TS) proce- 
dure to select only th e most v iable candidates for detailed 
analysis. In Paper I ([Kipping et al. 2012b). we discussed 
the three principal TS methods employed by HEK: i) au- 
tomatic target selection (TSA) ii) visual target selection 
(TSV) and iii) target selection opportunit ies (TSO). De- 
tails on all three methods are discussed in lKipping et al.l 
(|2012b( ). but this work will make use of TSA only. A 
dedicated TSV survey will be presented in a subsequent 
work. 

After the TSA stage, we also apply a target selection 
prioritization (TSP) selection process, which identifies 
the optimal targets for an exomoon hunt. 

2.1.2. Modifications to the TSA algorithm 

The automatic target selection (TSA) algorithm has 
been slightly modified since Paper I ([Kipping et al.l 
l2012bD . The main modification is to accommodate a 
continuous, and thus more realistic, minimum plane- 
tary mass estimation function. This mass function is 
required to estimate the maxi mum stable moon mas s 
around each KOI as described in lKipping et al.l ([2012b[ ). 
Previously, we considered three regimes: i) Super-Earths 
{Rp < 2.0 i?©) ii) Neptunes (2.0 i?e < Rp < 6.0 i?®) 
and iii) Jupiters {Rp > 6.0 i?©). The mass was esti- 
mated for Super-Earths using a terrestrial-scaling law 
from [Valencia et al.l (|2006j), whereas Neptunes were as- 
sumed to have a constant density of 1.7 g cm~^ for rea- 
sons discussed in iKipping et al.l ([2012bl ). Jupiters were 
not considered at all due to the hi gher potential for a 
false-positive (Sa nterne et al.l [20121) . These minimum 
mass estimates are required to evaluate the dynamical 
capacity of each KOI for hosting a moon and thus a min- 
imum mass provides a conservative lower limit. 

In this revised TSA algorithm, we make two major 
changes to the mass function: 1) we ensure a continuous 
mass-function 2) we allow this mass function to go into 
the Jupiter-regime. The first improvement is inspired 
by the fact that i n the Super-Earth/Neptune regime, 
the iValencia et all ([20061 ) mass function quickly exceeds 
10 M©, which leads to an abundance of TSA targets at 
this boundary (since higher mass planets have a better 
chance of hosting a moon). The second improvement al- 
lows us to consider the Jupiters as well and thus expand 
our se arch somewhat. 

The IValencia et all ([2006D expression of i? - AfO ^^ 
is invertible to M ~ R^-^, in units of Earth radii 
and masses. For R = 1.863 i?© the mass hits 10 M©, 
which we consider a sensible upper limit for Super-Earth 
masses. Recall that TSA is primarily interested in a con- 
servative estimate of the planetary mass in order to en- 
sure selected candidates have the best chance of being 
a good target. We therefore consider the Super-Earth 
regime to be modified by this new radius limit. To 
bridge the gap between Super-Earths and Neptunes in 
a continuous manner, we fix the mass to be 10 M© un- 
til (4/3)7rpNcptuno^P = 10 M©, which occurs for Rp = 
3.186 i?0 (we assume spherical planets). The region be- 
tween is dubbed "Mini-Neptunes" and such objects are 
assumed to have a mass of lOMfr,. 



Our previous boundary between Neptunes and 
Jupiters is also discontinuous. The mass of Jupiter- 
radius planets varies widely, but a sensible lower limit 
is to assume a Jovian bulk density. Much lower-density 
Jupiters do exist, so called inflated gas giants, but are 
thought to be due to their high irradiatio n environment 
since they are typically hot- Jupiters ([Burrows et al.l 
|2007[ ). TSA automatically excludes hot- Jupiters since 
they are too close to the ir star to maintain an exomoon 
([Weidner fc Horn^ 120101 ). We therefore adopt a Jovian 
density of 1.326 gcm~'^. To bridge the density discon- 
tinuity between Neptunes and Jupiters, we assume the 
bulk density linearly drops off between 6 i?© and 7 i?© 
to create a continuous mass function. Finally, for a Jo- 
vian density object, once the radius exceeds 10.963 i?© 
the mass will exceed a Jupiter masJ*^. For such cases, we 
set the upper limit on the mass to be 1 Mj. This is again 
in line with the conservative mass estimate requirements 
of TSA. In summary, we have: 

i) "Super-Earths"; < < 1.863 

ii) "Mini-Neptunes"; 1.863 < |j < 3.186 

iii) "Neptunes"; 3.186 < |^ < 6 

iv) "sub- Jupiters" ; 6 < -|^^ < 7 

v) "Jupiters"; 7 < |^ < 10.963 

vi) "super- Jupiters" ; 10.963 < ^ < oo 

These six regimes are described by the following mass- 
function: 



Mp 



where 
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if 3.186 < 1^ < 6, 
if 6 < < 7, 



if? < ^< 10.963, 
if 10.963 < -g^ < oo. 



PsubJup — PNcp 



1 - 



Rp-QRa 
Rm 



1 - 



PJup 
PNcp 



(1) 



2.1.3. TSA inputs 



The TSA algorithm can be executed for several dif- 
ferent inputs. Critically, one can choose the maximum 
orbital distance that an exomoon can reside at (in units 
of the Hill radius) /max, and the tidal dissipation fac- 
tor of the host planet Qp, both of which strongly affect 
th e maximum allowed exomoon mass via the expressions 
of iBarnes fc O'BrienI ([2002f l. TSAs are defined to sat- 
isfy the criterion that the host planet can maintain an 

This does not occur at 71, 492, 000 km= 1 Rj since this the 
equatorial radius of Jupiter and the planet is non-spherical. 
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Fig. 1. — Assumed mass function of Kepler Objects of Interest 
(KOIs) used in the updated Automatic Target Selection (TSA) al- 
gorithm of this work. The six regimes, from low-to-high radius, 
are i) Super-Earth ii) Mini-Neptune Hi) Neptune iv) sub-Jupiter 
v) Jupiter vi) Super- Jupiter. Expressions for mass function are 
provided in § \2.1.2\ The eight triangles represent the Solar Sys- 
tem planets. 



Earth-mass moon for 5Gyr (jKipping et al.ll2012b[) and 
thus these inputs affect the number of TSA candidates 
identified. 

An additional freedom is that one can alter the re- 
quirement of the signal-to-noise ratio (SNR) for a moon 
transit. We define SNR as an Earth-radius transit depth 
divided by th e combined differential p hotometric preci- 
sion (CDPP) ^ Christiansen et al.|[20l2[ ) over 6 hours (see 
Equation [2]) . Note that the CDPP values are taken from 
the B12 tables, as with all other TSA inputs. One can 
see that enforcing a higher SNR condition will naturally 
reduce the number of TSA candidates found. 



SNR- 



CDPPfi 



(2) 



Finally, one can choose whether we look at host planets 
in the Neptune-regime and smaller or whether we expand 
our search to include Jupiter-size d objects (which ten d 
to have a higher false-positive-rate. ISanterne et al.ll2012[ ). 
In this work, we treat quarters 1 to 9 from Kepler as the 
survey data used to look for exomoons. This leads to an- 
other criterion for TSA that Pp < 193.6 days, so that at 
least three transits exist in the Ql-9 Kepler photometry; 
a minimum requirement for exomoon searches. At the 
time of writing, QlO-13 have also recently become avail- 
able and we treat these data as follow-up photometry 
with application for particularly interesting candidates. 
Any candidate already selected as a target of opportu- 
nity (TSO) by the HEK project was not included in the 
TSA lists. We summarize the results of making these 
modifications in Table [TJ 

2.1.4. TSA results 

Curiously, Table[T]reveals that Jupiter-sized objects of- 
fer very little improvement in the number of viable candi- 
dates for moon hunting. This is certainly due to a dearth 
of Jupiter-sized objects in the Kepler-samjple more than 
anything else (B12). 

For regular satellites, such as the Galilean satellites 
around Jupiter, the formation of large moons is en- 



TABLE 1 

Summary of the number of candidates identified using TSA from 
the 2321 KOIs available for various inputs to the TSA algorithm. 
The left column modifies the maximum allowed exomoon distance 
from the planet, in units of the Hill radius, fma,x, for several 
reasonable guesses. The second column presents the results for 
two different assumptions of the tidal dissipation factor, Qp. The 
final column shows a vector-form of the number of candidates 
found for the assumption of {SNR > 1, SNR > 2, SNR > 3}. 



/max 


Qp 


Candidates found 


Ftp ^ 6 


0.9309 


lO'^ 


{274,92,34} 


0.4805 


10^ 


{91,35,13} 


0.3333 


10^ 


{39, 19,8} 


0.2500 


10^ 


{16,7,4} 


0.9309 


lO* 


{156,59,20} 


0.4805 


10-* 


{42, 20, 8} 


0.3333 


lO* 


{9,4,2} 


0.2500 


lO* 


{0,0,0} 


Rp > 6R® 


0.9309 


10^ 


{26,5,2} 


0.4805 


10^ 


{13,2,2} 


0.3333 


10^ 


{9,2,2} 


0.2500 


10^ 


{1,0,0} 


0.9309 


10"* 


{20,4,2} 


0.4805 


10-* 


{9,2,2} 


0.3333 


10-* 


{1,0,0} 


0.2500 


10-* 


{0,0,0} 



hanced by a massive primarv (ICanup fc Wardl 120061 : 
ISasaki et all [Ml IQgihara fc Idal 12012^. Hoover, i f 
the proposed mass scaling law of (|Canup fc Wardl[200l 
holds true, that Ms/Mp < 10~ ^, then such moons 
will be undetectable using Kepler (jKipping et al.l [2'0091. 
The moons we seek are therefore most likely irregu- 
lar satellites ar riving through capture or impact (e.g. 
Triton-Neptune 'Agnor fc H amiltonI I2006D . In such a 
case. Porter fc Gnmdy (2011^ argue that the mass of 
the primary has a much weaker effect with Neptunes and 
Jupiters retaining captured satelites with broadly equiv- 
alent efficiencies. It is therefore important to stress that 
Jupiter-sized KOIs do not hold a special significance for 
target selection over Neptunes. 

Due to the relatively small improvement in the overall 
number of candidates offered b y Jupiters, combined with 
their higher false-positive rate (jSanterne et al.l[2012t ) , we 
will not consider Jovian TSA candidates in this work, 
but we will return to them later in a future HEK survey. 

2.1.5. Selecting a TSA category 

Table [T] presents twenty- four different viable inputs for 
the TSA algorithm (for Neptunes or smaller). We must 
now select which input to use. The first point to bear 
in mind is that an excellent candidate will appear in 
multiple categories. For example, if a candidate satis- 
fies SNR> 3 then it will of course also satisfy SNR> 1. 
The task is therefore simply to move from the most con- 
servative estimate to the most optimistic and stop at the 
point at which we have the desired number of candidates. 

Due to computational constraints, we estimated we 
could analyze a handful of targets in this work. We there- 
fore attempted to select a category which yields around 
a dozen or so candidates and apply the final target se- 
lection prioritization (TSP) stage to filter out the best of 
those. 
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We choose to work with Qp = 10^ in what follows, 
since this option finds dramatically more high SNR sig- 
nals and thus the best chance for success. Next, we only 
consider candidates where the expected SNR> 2 to bal- 
ance between a high SNR and a significant number of 
candidates. We only consider KOIs with radii below 
Rp < 6i?0 (as reported by B12) for reasons discussed 
earlier. Finally, we opt for /max = 0.4805 which bounds 
all prograde satellites yet still returns 35 TSA candidates, 
which are listed in Table [5J 

2.2. Target Selection Prioritization (TSP) 
2.2.1. Overview 

The TSA algorithm has identified 35 KOIs as being 
suitable for an exomoon analysis. With this more man- 
ageable number, we can apply some more time intensive 
selection criteria as part of the target selection prioritiza- 
tion (TSP) process. In this work, we consider three TSP 
criteria, which sequentially increase in time requirements 
to evaluate: 

■ KOI must be in a single-transiting system 

■ SNR should hold-up when queried from the MAST 
archive 

■ KOI should not exhibit excessive time-correlated 
noise 

We discuss each of these three criteria in the following 
subsections. 

2.2.2. Multiplicity 

The first criterion eliminates KOIs which would require 
a more complicated and involved analysis due to their 
multiple nature. Multiples should induce transit timing 
variations (TTVs) on one another, which is also a sig- 
nature of exomoons. Eliminating these KOIs does not 
eliminate the possibility of planet-induced TTVs by any 
means (as r ecently demonstrated by the counter-example 
of KOI-872 lNesvornv et alll2012[ ). but it does make our 
task simpler. Later HEK surveys may relax this con- 
straint. Of the 35 TSAs, 11 were found to reside in mul- 
tiple transiting systems and were rejected. 

2.2.3. Cross-referencing CDPPs 

The TSA algorithm works by reading in a list of planet 
and star parameters for each KOI. The major source for 
such parameters comes from B12. This work also in- 
cludes estimates of the CDPP over 6 hours timescale, 
which TSA uses to estimate the SNR. During our inves- 
tigation, we noticed several cases where the CDPP values 
reported in the tables of B12 did not agree with those re- 
ported by MAST when queried. We therefore decided to 
cross-reference the SNRs calculated from the B12 CDPP 
to those given by MAST. 

In addition to the CDPP values differing between B12 
and MAST, KOIs with a low Tcff host star may yield 
unreliable Rj, es t imate s for reasons discussed in detail in 
iMuirhead et al.l (|2012[ ). These authors provide improved 
Ri, estimates for such systems, which we use where avail- 
able to compute a revised SNR. Since MAST provides 
the CDPP values for each quarter, we evaluate the mean 



and standard deviation of the SNR across all LC quar- 
ters. 

We find that 12 of the remaining 24 KOIs have a mean 
SNR< 2.0 when we used the revised values and these 
candidates are summarily rejected. This leaves us with 
12 KOIs. 

2.2.4. Removing KOIs with excessive time- correlated noise 

With 12 KOIs remaining, we are now ready to con- 
sider the most time-consuming TSP test. In the limit of 
a perfectly well-behaved star and instrument, the noise 
should be purely due to photon noise and thus behave 
as a Poisson distribution. Since the number of photons 
is large, the noise is very well described as a Gaussian 
distribution and has no frequency dependency; so-called 
"white noise" . 

Time-correlated noise refers to noise which has the 
property that the probability distribution of values for a 
given measurement is not independent of previous mea- 
surements. This is a problem for the HEK project since 
time-correlated noise can mimic dips, bumps and dis- 
tortions due to an exomoon. Whilst many methods ex- 
ist to tackle time-correlated noise, they require various 
assumptions about the data's behavior and invariably 
greater computational overhead. Since a significant frac- 
tion of Kepler 's targets have their photome try dominated 
by uncorrelated noise (| Jenkins et al.|[20Tol ). the simplest 
strategy to deal with time-correlated noise is reject any 
KOIs exhibiting an excess on the timescale of interest. 
On timescales of days to weeks, one invariably finds time- 
correlated fiux modulations, which could be considered a 
form of time-correlated noise, typically due to focus drift 
or stellar rotation. However, the timescale of interest for 
transiting exomoons is of order-of-magnitude one hour. 
Therefore, these long term variations do not affect our 
analysis and should be detrended out appropriately (see 

At this timescale of interest, it makes no difference to 
us whether excessive time-correlated noise is of instru- 
mental or stellar origin since we have no intention of at- 
tempting to correct for it. Instead, our strategy is simply 
to reject all candidates showing excessive time-correlated 
noise. The question then becomes, what do we define as 
excessive time-correlated noise? 

There is a dizzying number of metrics at our dis- 
posal for this task and we here seek a simple, com- 
putational efficien t expre ssion. A classic metric is the 
IDurbin fc Watsoiil ()1950l ) statistic, d, which uses auto- 
correlation to test whether a time series is positively or 
negatively autocorrelated. The Durbin- Watson statistic 
is given by 



where are the residuals and N is the number of data 
points. The value of d always lies between and 4, with 2 
representing an absence of autocorrelation, d < 2 repre- 
senting positively-autocorrelated noise (expected for in- 
strumental/astrophysical sources) and d > 2 represent- 
ing negatively-autocorrelated noise (anomalous and un- 
physical; we do not expect to see a significant excess of 
this). 
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TABLE 2 

Kepler candidates identified for an exomoon search using TSA. SNR is defined in Equation\^ The multiplicity denotes the number of 
KOIs in each system. Bold and bold-italic highlighted rows are those targets chosen for target selection prioritization (TSP) with the 
former being those which we accepted by TSP and the latter being those rejected, f = KOI dropped by Kepler-team. * = Durbin- Watson 
statistic (EauationWi indicated high (> 3cr) probability of correlated noise at 30 minute timescale in both PA and PDC data and the 

candidate was consequently rejected. 



KIC 


KOI 


Pp [days] 


SNR (B12) 


SNR (MAST) 


Multip 


3425851.01* 


268.01 


110.4 


5.76 


6.65±0.36 


1 


11623629.01 


365.01 


81.7 


6.39 


5.45 ± 0.62 


1 


7296438. 01 


364.01 


173.9 


4.68 


4.57±0.19 


1 


5966322.01 


303.01 


60.9 


2.62 


3.18 ±0.35 


1 


8292840.02 


260.02 


100.3 


2.69 


3.29 ±0.27 


2 


9451706.01 


271.01 


48.6 


2.32 


2.96 ±0.38 


2 


9414417.01* 


974.01 


53.5 


2.28 


2.87±0.25 


1 


9002278.03 


701.03 


122.4 


2.18 


2.77 ±0.11 


3 


9965439.01 


722.01 


46.4 


2.22 


2.31 ± 0.19 


1 


11622600.01 


1876.01 


82.5 


3.17 


2.27± 0.10 


1 


7199397.01* 


75.01 


105.9 


2.61 


2.22±0.11 


1 


11297236.01 


1857.01 


88.6 


2.10 


2.19 ±0.09 


1 


9349482.01* 


2020.01 


111.0 


3.72 


2.17±0.22 


1 


10810838.01 


174.01 


56.4 


2.04 


2.15 ± 0.16 


1 


10471621.01 


2554.01 


39.8 


3.19 


2.10 ±0.33 


2 


7761545.01 


1472.01 


85.4 


2.60 


2.04 ±0.35 


1 


8686097.01 


374.01 


172.7 


2.09 


1.92 ±0.07 


1 


12121570.01 


2290.01 


91.5 


2.73 


1.88 ±0.07 


1 


9661979.01 


2132.01 


69.9 


2.11 


1.87 ±0.14 


1 


2449431.01 


2009.01 


86.7 


3.82 


1.81 ±0.25 


1 


2443393.01 


2603.01 


73.7 


2.46 


1.79 ±0.07 


1 


5526717.01 


1677.01 


52.1 


2.65 


1.77 ±0.16 


2 


10027247.01 


2418.01 


86.8 


3.66 


1.70 ±0.08 


1 


11037335.01 


1435.01 


40.7 


2.22 


1.54 ±0.07 


2 


12400538.01 


1503.01 


76.1 


2.00 


1.20 ±0.22 


1 


10015937.01 


1720.01 


59.7 


3.98 


1.07 ±0.07 


1 


11656918.01 


1945.01 


82.5 


2.31 


1.06 ±0.08 


2 


11176127.03 


1430.03 


77.5 


3.18 


1.05 ±0.10 


3 


8611781.01 


2185.01 


77.0 


2.34 


1.05 ±0.11 


1 


8892157.02 


2224.02 


86.1 


2.92 


1.01 ±0.08 


2 


6765135.01 


2592.01 


175.6 


2.53 


0.80 ±0.05 


1 


8758204.01 


2841.01 


159.4 


3.49 


0.75 ±0.08 


1 


9030537.01 


1892.01 


62.6 


4.18 


0.71 ±0.10 


1 


8240904.02 


1070.02 


107.7 


2.66 


0.40 ±0.09 


3 


8240904.03 


1070.03 


92.8 


2.66 


0.40 ±0.09 
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In calculating d, there exists a degree of freedom re- 
garding over what cadence should we evaluate the statis- 
tic i.e. what timescale do we consider most relevant for 
an exomoon search? An exomoon transit or distortion 
could occur on a timescale of a few minutes or a few hours 
and so we selected 30 minutes for the simplicity that no 
binning is required for LC data and that the timescale is 
consistent with exomoon features. The next question is 
what value of d should one consider acceptable? 

Although test statistics are available via lookup tables 
for instance, these statistics assume regularly spaced time 
series which we often do not have for Kepler data, mostly 
due to outlier rejections and data gaps. Instead of using 
such statistics, we generate 1000 Monte Carlo simula- 
tions of Gaussian noise for the exact time sampling of a 
given data set to reproduce the expected posterior distri- 
bution of d for data with no time-correlated noise. The 
Gaussian noise is generated assuming each data point is 
described by a Gaussian distribution with a standard de- 
viation given by its associated uncertainty and a mean 
of unity. 

The quantity d is only calculated on data locally sur- 
rounding transits to within twice the timescale of the Hill 



sphere, 2rHiu (with the planetary transit itself excluded), 
since only this data is relevant for our exomoon hunt. We 
compute d for each and every transit event and then take 
the mean over all transit epochs, d. The same process is 
applied to the 1000 synthetic time series, which behave 
as Gaussian noise with the exact same cadence and time 
sampling as the original data. We also apply the same 
final-stage local linear detrending used on the real data 
to every synthetic data set (as is described later in @. 
The 1000 synthetic time series are converted into 1000 d 
metrics in the same as the original data and this is used 
to compute a probability distribution of d in the case of 
Gaussian noise. 

We then compare the real d metric with the simulated 
distribution to evaluate whether our data set is consis- 
tent with a lack of autocorrelated noise. Any KOIs which 
show > 3 (7 autocorrelation (as determined by the d met- 
ric) in both the PA (Photometric Analysis) and PDC 
(Pre-search Data Conditioning) detrended LC data are 
rejected (SC data is not used in this selection phase but 
is used later). We do not anticipate that this will re- 
move genuine moon signals since such events would be 
temporally localized rather whereas autocorrelation at a 
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TABLE 3 

Durbin- Watson (d) statistics of TSAs. Each vector displays two 

numbers; the first for the PA data and the second for the 
PDC-MAP. Kepler data usually yields d < 2 indicating positive 
serial correlation, as would be expected for some hidden 
systematic error source. Rows in italics are those rejected for 
having excessive autocorrelation. 



KOI 


d of LC 


(1-FAP) of autocorrelation (cr) 


KOI- 364. 01 






KOI-303.01 


{1.922, 1.893} 


{2.0,2.6} 


KOI- 974. 01 


{1.423,1.498} 


{11.9,10.8} 


KOI-268.01 


{1.446,1.415} 


{9.5,10.5} 


KOI-1472.01 


{1.951, 1.945} 


{1.4,1.4} 


KOI-722.01 


{1.951, 1.944} 


{1.7,1.9} 


KOI-365.01 


{1.905, 1.847} 


{2.2,3.1} 


KOI-174.01 


{1.947, 1.999} 


{1.5,0.6} 


KO 1-75. 01 


{1.422,1.471} 


{13.6,13.4} 


KOI- 2020. 01 


{1.643,1.633} 


{6.4,6.7} 


KOI-1857.01 


{1.901, 1.902} 


{2.3,2.3} 


KOI-1876.01 


{1.909, 1.930} 


{0.5,0.4} 



30 minute timescale must be present throughout the en- 
tire time series (for the particular transit epoch under 
analysis). After applying this test to each target, we find 
that only 7 of the 12 remaining KOIs pass this test, as 
listed in Table H Note that all 12 KOIs were fully de- 
trended in exactly the same way, as is described in fj3l 
However, only 7 of these 12 are actually fitted with a 
transit light curve model - the most resource intensive 
stage of the entire process. 

3. DETRENDING THE DATA WITH CoFiAM 

Data is detrended using a custom algorithm which we 
dub Cosine Filtering with Autocorrelation Minimization 
(CoFiAM), which is described in this section. 

3.1. P re- Detr ending Cleaning 

In all cases, we performed the detrending procedure 
twice; once for the PA data and once for the PDC-MAP 
data. In what follows, each transit is always analyzed 
independently of the others i.e. we obtain a detrended 
light curve unique to each transit event, not each quar- 
ter. The first step is to visually inspect each quarter and 
remove any exponential ramps, flare-like behaviours and 
instrumental discontinuities in the data. We make no at- 
tempt to correct these artefacts and simply exclude them 
from the photometry manually. We then remove all tran- 
sits using the B12 ephemerides and clean the data of 3 cr 
outliers from a moving median smoothing curve with a 
20-point window (for both LC and SC data). 

3.2. Cosine Filtering with Linear Minimization 
(CoFiAMj 

The remaining unevenly spaced data is then regressed 
using a discrete series of harmonic cosi ne functions, 
which act as a high-pass, low-cut filter ()Ahmed et al.l 
119741) . The functional form is given by 



N„ 



fk{ti) = oo + ^ 



Xfe sm 



2'Ktik 
2D 



Vk cos 



/■2'!rtik\ 
\ 2D ) 



(4) 



where D is the total baseline of the data under anal- 
ysis, ti are the time stamps of the data, Xk & Uk are 



model variables and iVordor is the highest harmonic or- 
der. Equation |4] may be more compactly expressed as 
a cosine function with a phase term, but the above for- 
mat illustrates how the equation is linear with respect 
to Xk and t/fe. This means that we can employ weighted 
linear minimization, which is not only computationally 
quicker than non-linear methods, but also guaranteed to 
reach the global minimum. In our regression, the data 
are weighted by the inverse of their reported standard 
photometric errors. 

Harmonic filtering has b een previously used to correct 
CoR oT Ce.g. iMazeh et a l. 2010) and Kepler light curves 
fe.g. iKipping fc Bakosll2011a: Kipping fc Bakosll2011bl: 
IKipping fc Spiege]|l201lD and is attractive for its simplic- 
ity, computational efficiency and ability to preserve the 
transit shape. 

3.3. Frequency Protection 

There are many possible choices of A^ordcr, but above 
a certain threshold the harmonics will start to appear at 
the same timescale as the transit shape and thus distort 
the profile, which is undesirable. The transit light curve 
of a planet can be considered to be a trapezoid to an ex- 
cellent approximation, fo r which an analyt ic Fourier de- 
composition is available. iWaldmamil (|2012f ) showed that 
a equatorial trapezoidal transit light curve is described 
by the following Fourier series: 



/(t) = 8V2(5 sin(l/Ti4) + 



sin(3/Ti4) sin(5/ri4) 
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(5) 



Under the approximation of a trapezoidal light curve 
then, the lowest frequency is thus (l/ri4). Another way 
of putting this is that the highest periodicity is ri4 (i.e. 
the transit duration) and so if we protect this timescale, 
and all shorter timescales, the transit light curve should 
be minimally distorted. Let us therefore choose to pro- 
tect a timescale xTn, where x is a real number greater 
than unity and T14 is the first-to-fourth contact duration 
reported in B12. The timescale xTn can be protected 
by imposing 



2D 

AxTu 



(6) 



Ideally, one would wish to impose a; 3> 1 in all cases 
to provide some cushion, but in reality such a condi- 
tion means A^ordor is small and the ability of the regres- 
sion algorithm to obtain a reasonable fit to the data be- 
comes poor. In contrast, going to higher values of A'ordor 
leads to a better regression in the sense, but increases 
the risk of higher harmonics distorting the transit pro- 
file since strictly speaking we require x 3> 1. Therefore, 
there exists a trade-off between these two effects and one 
might expect an optimal choice of x to exist for any given 
data set. Selecting such an optimum requires a quanti- 
tative metric which we aim to optimize. 

3.4. Autocorrelation Minimization 

For any optimization problem, one must first define 
what it is we wish to optimize or minimize. In this work, 
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we identify the primary objective to be that the transit 
hght curve contains the lowest possible degree of time- 
correlated noise around each transit, which could lead to 
false-positive moon signals. We therefore require some 
metric to quantify the amount of autocorrelation and op- 
timize against. As discussed earlier, the Durbin- Watson 
statistic is a useful tool to this end and evaluates the de- 
gree of first-order autocorrelation in a time series. Our 
objective is therefore to choose a value of x such that 
the Durbin- Watson statistic is consistent with the low- 
est quantity of autocorrelation, when evaluated on the 
data surrounding the planetary transit. Another way of 
putting this is that we wish to choose the value of iVordcr 
which minimizes {d — 2)^, when evaluated on the data 
within 2THiii of the time of transit minimum (where THiii 
is the Hill timescale). 

Before we can begin our optimization search, one must 
define the allowed range of A^ordcr through which we can 
search. Recall that our expressions protect a timescale 
xTi4 and thus one must choose a; > 1 in order to not dis- 
turb the transit profile. Consequently, we chose the low- 
est allowed value of x to be 3, such the timescale 3Ti4 is 
never perturbed. This factor of three cushion is to allow 
for transit duration changes, leakage of the harmonics 
for real transits and longer exomoon transits. Protecting 
this timescale corresponds to a maximum allowed value 

of iVordcr of 



border, max -, ^-.rj^ • [J ) 

Our detrending algorithm, CoFiAM, regresses the Ke- 
pler time series to the harmonic series given by Equa- 
tion 2] in a least squares sense and repeats this regression 
for every possible integer choice of A'ordcr between 1 and 
A'ordor.max- In this way, wc cxplore dozens of different 
regressions which all satisfy the conditions of providing 
a good least squares fit to the light curve and protect a 
timescale > 3Tn. We then simply scan through the final 
list of d values and define the optimal detcnding func- 
tion to be the harmonic order which minimized (d — 2)^ 
when evaluated on the out-of-transit data within 2THin of 
the eclipse. If a quarter contains more than one transit, 
we always repeat the entire procedure for each transit 
to ensure the data associated with the transit is fully 
optimized for our exomoon hunt. 

As before, the d statistic is computed on a timescale of 
30 minutes which means that no binning is required for 
the LC data. For the SC data, we bin the data up the 
long-cadence data rate. 

Once the optimal detrending function has been found, 
we divide all data within 2THiu (including the transit) by 
fk{tk) to correct for the long-term variations. We also 
apply a second outlier rejection of 10 a filtering (to allow 
for unusual anomalies in the transit) from a moving 5- 
point median. For short-cadence data, we instead use a 
3 a filtering on a 20-point moving median. In some cases, 
we relaxed this outlier rejection when we felt the filter 
was removing potential exomoon signals. We clip out 
the data within ±2THiii and apply a linear fit through the 
out-of-transit data to remove any residual trend, which 
acts as a final normalization. The process is repeated 
for all transits and the surviving light curves are stitched 
together to form a single input file for our fitting code. 



Any transit epoch with > data points in the window 
r ± 2THiii is accepted into the final file. Additionally, the 
final file always uses SC data over LC data, where such 
data exists. 

4. LIGHT CURVE FITS 
4.1. Overview 

Model light curves of a transiting pl anet are g e nerate d 
using the LUNA algorithm described in iKippinei ()2011af) . 
LUNA is an analytic photodynamic light curve modeling 
algorithm, optimized for a planet with a satellite. LUNA 
accounts for auxiliary transits, mutual events, non-linear 
limb darkening and the dynamical motion of the planet 
and its satellite with respect to the host star. In the 
case of a zero-radius and zero-mass sa tellite, the LUNA 
expressions are equivalent to the familiar [Mandel &: Agoll 
(2002) algorithm. 

For any given model, we regress the data to the model 
parameters using the MultiNest algorithm (Fero z et al.l 
I2009a!b[). MultiN est is a multimodal nested sampling 
(see ISkillindl200l algorithm designed to calculate the 
Bayesian evidence of each model regressed, along with 
the parameter posteriors. By comparing the Bayesian 
evidences of different models, one may conduct Bayesian 
model selection, which has the advantage of featuring 
a built-in Occam's razor. For each KOI in our survey, 
we always regress the following models as a minimum 
requirement: 

■ Planet-only model with variable baselines using 
theoretical quadratic limb darkening coefficients; 
model Vp. 

■ Planet-only model with variable baselines and free 
quadratic limb darkening coefficients; model Vp^ld- 

■ If logZ(Vp,LD) > log2^(Vp), then we repeat Vp fix- 
ing the LD coefficients to the maximum a-posteriori 
LD coefficients from model Vp,ld (the best LD 
coefficients will be adopted throughout from this 
point on), in a fit which we dub as Vp^map- 

■ Planet-only model with fiat baseline over all 
epochs, J^p. If logZ(J"p) > logZ(Vp) then we 
use a flat baseline in all following fits (to reduce 
the number of free parameters), which was found 
to be always true thanks to CoFiAM employing a 
final-stage normalization. 

■ Planet-only model with variables times of transit 
minimum (i.e. TTVs); model J^ttv- 

■ Planet-only model with each transit possessing 
unique transit parameters to allow for both TTVs 
and TDVs; model Vy. 

■ Planet-with-satcllitc fit; model Ts- 

m Planet-with-satellite fit assuming a zero-mass 
moon; model J^s,MO- 

■ Planet-with-satellite fit assuming a zero-radius 
moon; model J^s.RO- 
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The question as to why we switch from local baselines 
to a global baseline is discussed in the next subsection, 
In all fits, we use the same data set throughout, 
which is usually the PA data. This is because the PDC- 
MAP data is subject to numerous detrending processes 
that do not necessarily preserve exomoon signals. How- 
ever, if the PA data yields a d statistic with more than 
3(7 confidence of autocorrelation but the PDC-MAP is 
below 3cr, then the PDC-MAP data is used in the fits 
instead. However, Table [3] reveals how there is no such 
instance in the sample of KOIs studied in this work. 

4.2. Planet-only Fits 

The first stage of our fitting process always begins with 
planet-only fits. The purpose of these fits is to i) verify or 
obtain reliable limb darkening parameters ii) serve as a 
baseline for comparison with the planet-with-moon fits. 
Initially, we employ fixed limb dark ening co e fficien ts, cal- 
culating theoretical values from a iKurucd (|2006[ ) style- 
atmosphere integrated over the ifep/er-bandpass. The 
computation is performed by a code written by I. Ribas 
and associated details can be found in lKipping fc BakosI 
((2011a). 

In these fits, the five basic transit parameters are p, 

b, pV^ , P and r. The choice of these five parameters 
is fairly commonplace in the exoplanet literature, except 

2/3 

for perhaps p+ . This parameter is used so that the pos- 
teriors of pV^ have a uniform prior and can be utilized 
with the Multibody A sterodensity Pr o filing ( MAP) tech- 
nique, as discussed in IKipping et al.l (|2012al ). Although 
it is not the purpose of this work to conduct MAP, the 
posteriors are available upon request so that these stud- 
ies can be facilitated in the future without re-executing 
the light curve fits. 

For all KOIs, an estimate of P and r exists from B12 
which we use to define a uniform prior of ±lday either 
side of the estimate for both P and r. The other three ba- 
sic parameters also have uniform priors ofO<p< 1,0< 
5 < 2 and 7.6499 kg^/^ m-^ < p2/3 ^ 6097.85 kg^/^ 
(covering the main-sequence of stars between spectral 
types of M5 to FO with a factor of 10 cushion at each 
boundarv lCox|[2000l ). The transit epoch, r, can be cen- 
tered on any one of the transits in the Ql-9 time series. 
We choose an epoch which is nearest to the median time 
stamp of Ql-9 data, in order to minimize degeneracy be- 
tween P and T in the fits (jPal 2009). 

Whether the fit is using a variable baseline (Vp) or 
a flat baseline (J-p), we adopt a uniform prior of 0.95 < 
OOT < 1.05. For F-p model, the total number of free pa- 
rameters to be marginalized over (and thus explored by 
MultiNest) is just 5+1=6. The Vp model has 5+iV free 
parameters. Experience with MultiNest shows that fit- 
ting models with more than >20 free parameters becomes 
dramatically more time-consuming and so the flat base- 
line model is useful later for exomoon fits with a greater 
number of basic parameters. It is for this reason that we 
transition from local baselines to a global baseline as the 
model complexity increases. 

When fitting for limb darkening parameters, we use 
quadratic limb darkening so that only two degrees of free- 
dom are required, yet the c urvature of t he light curve 
can be modelled effectively (jClaretl |2000[ ) . For consis- 
tency, we employ quadratic limb darkening when we use 



the fixed limb darkening parameters too. We fit for the 
terms ui and ui-\-U2 since they are bounded by the physi- 
cally motivated low er and up per limits of < mi < 2 and 
< Ui -\-U2 < 1 (see lCarter e t al. 2009' and lKipping et al.l 
12012b'). 

In both the planet-only fits, and all subsequent fits, 
we account for the integration tim e of the long-ca dence 
data using the resamphng method (lKippinell2010aD with 
jVrBsa.Tp = 5 (as generally recommended in iKippind 
I2010al) . AU MultiNest fits wi ll also e i nploy 4 000 live 
points, as recommended by iFeroz et al.l (|2009bl ) for evi- 
dence calculations. 

4.3. TTV & TDV Fits 

A transit timing variation (TTV) fit is performed by 
assigning each transit epoch a unique time of transit min- 
imum, Ti. All other parameters are kept global as before 
in J^p. This TTV fit is dubbed J^-tv- The period prior 
is changed to a Gaussian prior assigned from the poste- 
rior of P from the model J^p fits. Without a constraining 
prior on P, allowing every transit epoch to have a unique 
transit time would mean that P would be fully degener- 
ate with the Ti parameters. 

Transit duration variations (TDVs) may be due to ve- 
locity changes or impact-parameter changes of the ob- 
served planet. These changes induce not only changes 
in the transit duration, but changes in the derived a/i?+ 

(and thus pV^) and b. We also search for changes in the 
apparent p value due to spot activity or exomoon mutual 
events, for example. Additionally, TDVs are expected to 
occur in dynamic systems exhibiting TTVs too. With 
so many degrees of freedom required, the fits would cer- 
tainly involve a large number of free parameters making 
the regression very time consuming with MultiNest. 
To solve this, one may paradoxically increase the num- 
ber of degrees of freedom again by allowing for variable 
baselines (OOT). By doing so, all six transit parame- 
ters are independent for each epoch (i.e. there are no 
fitted global parameters) and thus the fits can be con- 
ducted on each epoch separately and then the sum of the 
log Bayesian evidences will give the global log Bayesian 
evidence. These individual transit fits are very fast to 
execute and may be run simultaneously. For these rea- 
sons, we dub the regression Vy meaning variable transit 
parameters for each epoch and the subscript V denotes 
variable baselines too. 

TDVs may be defined in several ways, unlike TTVs 
which have less ambiguous definition. In this paper, we 
define the TDVs to be the variation of the parameter 
rB*/2, where Tb* is the duration for the planet's cen- 
ter to enter to stellar disc and subsequently leave. We 
use this definition rather than the first-to-fourth contact 
duration, for example, since Tb* has the lowest rela- 
tive u ncertainty wh en the limb darkening coefficients are 
fixed ([Carter et al.l 120081. We divide this duration by 
two since the theoreti cal uncertainties o n Tb* are ex- 
actly twice that of t ([Carter et al.l 120081 ) and thus our 
derived TTVs and TDVs should exhibit similar scatter 
and scale, which makes for useful comparisons. 

4.4. Planet-with-Moon Fits 

In general, we make the following assumption: ex- 
omoons may be randomly oriented but have nearly 
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circular orbits due to tidal dissipation (e .g. see 
IPorter fc Grll^lMl iHeller fc Earned [20ll . There- 
fore, our survey-mode fits do not consider fitting for the 
orbital eccentricity of the exomoon. We find that this 
dramatically improves the speed and stability of our fits 
using MultiNest. 

In contrast to exomoons, there is no reason to expect 
exoplanets to have zero eccentricity, especially at long pe- 
riods, due to the much longer circularization timescales 
and the possibility for planet-planet forcing or Kozai mi- 
gration. Despite this, we choose to assume a circular 
orbit in the survey light curve fits in this paper. The 
advantage of doing so is firstly to save MultiNest ex- 
ploring an additional two free parameters and secondly 
to save solving Kepler's equation numerically at every 
time stamp. 

We justify our choice on the basis that the maximum 
stable orbital separation of a moon decreases rapidly with 
res pect to the host planet' s orbital eccentricity, as shown 
bv iDomingos et aTT ()20Q6[ ). Here, the authors find the 
maximum separation scales as ^ (1 — ep). Further, the 
maximum stable exomoon mass around a host planet 
sc ales as this separation to the index of 13/2, as shown 
bv lBarnes fc O'BrienI (jl002). We therefore expect that 
the maximum moon mass around a host planet scales as 
^ (1 — ep)^'^/^. On this basis, an eccentricity of even 0.1 
halves the maximum stable exomoon mass and an eccen- 
tricity of 0.3 reduces it by an order-of-magnitude. Future 
HEK surveys may explore eccentric planet solutions, but 
for this paper computational constraints limit our survey 
to circular systems for the reasons discussed. 

Another assumption we make is that exomoons orbit in 
a prograde sense. The gravitational influence of a satel- 
lite induces transit timing variations (TTV) and velocity- 
induce d transit duration variations (TDV-V) ([Kippinel 
l2009a| ). Both of these effects are insensitive to the sense 
of the moon's orbital motion. Additionally, satellites in- 
duce transit impact parameter induced transit duration 
variations (TDV-TIP), which typically have an ampli- 
tude of around an order-of-magnitude less than that of 
TDV-V. Howev er, TDV-TIP is sensitive to the sense of 
orbital motion (jKippingl l2009bl ). Therefore, by treat- 
ing all exomoons as prograde in our survey, a retrograde 
moon would have not its TDV-TIP effect modeled cor- 
rectly and thus an implicit assumption is therefore that 
TDV-V>TDV-TIP. As with the eccentricity assump- 
tion, this allows us to halve the parameter volume to be 
scanned through and thus expedite the fitting procedure. 

In addition to the standard planet-with-moon model, 
we try two "unphysical" fits where we fix the moon's 
mass and then the radius to be zero. These fits are use- 
ful in the vetting stage since a moon detection should 
not yield an improved Bayesian evidence with unphysi- 
cal properties, such as zero-radius. These fits are dubbed 
•^S,Mo and J-^s,Ro for the zero-mass and radius cases re- 
spectively. 

A careful choice of the parameter set and priors is cru- 
cial for the moon's parameters, since the expected signal- 
to-noise is low and thus priors can be expected to play 
an increasingly signific ant role in the derived results. In 
IKipping et alT()2012bl ). we suggested Ms/Mp, Rs/Rp, 

Psj Pp , 4>S, *S and fts with uniform priors for all. In 
this work, we have found that these priors were not fully 



adequate and our greater experience has led us to pro- 
pose a modification. Firstly, Pg is now fitted with a Jef- 
frey's prior since it spans several orders of magnitude and 
the low periods require dense sampling due to the bunch- 
ing up of harmonics. Secondly, we have exchanged is for 
cos is to impose an isotropic prior. Thirdly, we have ex- 
changed p^p^ for asp/Rp (the separation between the 
planet and moon in units of the planetary radius). This 
last change is geometrically motivated and means that 
MultiNest scans for moon transits evenly in time and 
space from the primary planet event. Except for Ps, all 
terms have uniform priors. 

The parameter cosig has the intuitively obvious prior 
of U{—1, 1}. Similarly, we use U{—tt,'k} for (jis and Vis 
with the exception that the boundary conditions are pe- 
riodic and thus the parameters are considered "wrap- 
around" parameters in MultiNest. The mass and ra- 
dius ratios have the very simple boundary conditions of 
being between zero and unity i.e. U{Q, 1}. 

The upper boundary condition on Ps is given by 
Pp/V^, w hich repr e sents t he edge of the Hill sphere as 
proved by iKippingI (|2009aD (i.e. /n,ax < !)■ T h e lower 
boundary is less obvious and in IKipping et al.l (|2012bf ) 
we proposed ^ 2 hours as a rough estimate. Since 
< Rs / Rp < 1, the maximum size of the moon is Rp. 
If this is the case, then the closest separation allowed 
before contact would occur is 2Rp. On this assumption, 
o ne may derive a lower limit for Ps- From Equation 7 
of IKippind (|2010bl ). one can re-arrange to make Ps the 
subject: 

_ (37r)i/2(a5g/P,)3/2(l + Mg/Mp)(Mg/Mp)i/^ 

(8) 

Substituting ass with a5p/(l -I- Ms/Mp) and then 
replacing asp with the minimum allowed value of 2i?p 
we have: 

P - /24^ f P ^ 3/2 (1 + Ms/Mp){Ms/MpY/^ 

(9) 

Exploiting the fact ps — pp{Ms/Mp){Rs / Rp)^^ and 
cleaning up the expression we find: 



/24^ / 1 

Ps,„.in - V y pp{l + Ms/Mp) ^^^> 

The maximum allowed value of (1-1- Ms/Mp) is 
2 and the maxi mum physically plausib le value oi pp 
is 27950 kg m-3 (|Kipping et ai.l l2012b[ ). This yields 
Ps.min = 0.0520311 days (1.25 hours) which we employ 
as our lower boundary condition for Ps- 

For asp/Rp, the lower boundary condition is simply 
2, which for < i?s/Pp < 1 guarantees no contact (for 
a circular satellite orbit). The maximum requires an- 
other small derivation. We consider the maximum to be 
the Hill sphere, i?Hiii — o-b* 

[Afp/(3Af*)]i/3; therefore we 

have: 
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/asp\ _ {aB*/R*) / A/p \ 1/3 

Vi?p /max ^ p 

Replacing the mass terms with densities, we obtain 



(11) 



/ asp 
V Rp 



0-1/3 1/3 -1/3/ /n X 

= 3 ' ppp p* [aB*/R*)- 



(12) 



If we assume {Mp + Ms) <C Af*, then = 
{3n{aB* / R*)^) / {GPg) which can be substituted into the 
above expression to give: 



asp \ 
Rp /, 



32/37^1/3 



,2/3 



(13) 



To estimate this value, we adopt the maximum al- 
lowed planetary den sity of pp^max = 27950 kg m"'^ 
(|Kipping et al.ll2012b[ ) and use the maximum a-posteriori 
value of Pb from the planet-only fits (technically model 
•^p)- OLSP I Rp is then fitted with a uniform prior between 
2 and this maximum value. 

Finally, we instruct MultiNest to ignore any trials 
which yield an unphysical density for the planet, which 
we consider to be bet ween 80kgm~'^ and 27950 kg m^'^ 
(|Kipping et al.ll2012b( l. We apply the same constraint to 
the satellite except the lower allowed limit is Okgm^"^. 
This is imposed so that zero-mass moons can be explored 
meaning that in the case of a null-detection the posterior 
of Ms/Mp can still reach zero. For models, J^s.mo and 
•^S.RO we deactivate the constraint on the satellite density 
since the satellite is specifically defined to be unphysical. 

4.5. Detection Criteria 

An exomoon has never been detected and thus one is 
forced to seriously consider what constitutes a "detec- 
tion" in such a new area. The signal-to-noise will in- 
evitably be at, or close to, the limit of Kepler, the signal 
will vary in phase and time and may manifest simply as 
a slight distortion to a planet's transit profile. If one 
fits a planet-with-moon model to real data, the extra de- 
grees of freedom will inevitably lead to an improved 
relative to a planet-only fit. Clearly an improved is 
not sufficient to claim an exomoon has been found. This 
concern was one of the driving reasons why the HEK 
project adopted Bayesian model selection (Kippin g et al.l 
I2012b( l available through MultiNest, since such com- 
parisons implicitly penalize models for using extra pa- 
rameters (Occam's razor). However, even the Bayesian 
evidence is not a tool which can be wielded blindly to 
claim exomoon detections. 

Although our detrending process CoFiAM minimizes the 
amount of autocorrelation, the data will always possess 
some quantity of time-correlated noise. The likelihood 
function employed by LUNA is a Gaussian likelih ood ex- 
pression (see Equation 20 of lKipping et al.ll2012bD and so 
this assumption will never be strictly true. It is therefore 
possible, and in fact quite common, that the Bayesian 
evidence of a planet-with-moon fit will be superior to a 
planet-only fit for an isolated planet, even with the built- 
in Occam's razor of Bayesian model selection. A pos- 
sible remedy would be to employ a more sophisticated 
likelihood function but the computational demands of 



MultiNest make this unrealistic as fits typically take 
weeks to run even on modern clusters. Instead, we stress 
that a superior Bayesian evidence is not tantamount to a 
detection with exomoon fits since LUNA can generate dis- 
tortions both in- and out-of-transit which can describe 
certain time-correlated noise features. 

We therefore consider that a superior Bayesian evi- 
dence of a planet-with-moon fit to a planet-only fit (i.e. 
logZ(J's) > log2(J"p)) is a requirement for a "detec- 
tion", but no t a pro of in of itself. As discussed in 
iKipping et al.l ()2012b( ). we set the significance level of 
this improvement to be in 4 tr or greater in order to qual- 
ify. This discussion therefore indicates that other detec- 
tion criteria are required. 

One of the easiest tests is that the posteriors s hould be 
physica lly plausible. For example, as shown in iKippinj 
(l2010bD . exomoons allow us to measure the ratio Afp/Af* 
and which for an assumed Mi, yields A^fp. It is easy to 
check whether the derived A^p is consistent with the de- 
rived Rp from known planet populations. For example, 
a candidate planet yielding a Jupiter-mass and an Earth- 
radius can be easily dismissed. 

Further, a planet-with-moon model must be superior 
to both unphysical moon models considered i.e. the zero- 
mass moon and the zero-radius moon models. If the 
zero-mass moon model is superior, one should suspect 
starspots or correlated noise to be responsible. If the 
zero-radius moon is superior, one should suspect TTVs 
from a non-moon origin. In both cases, it is also possible 
that the signal-to-noise of the moon signal is presently 
too low to make a confirmed detection (but more data 
may change this). 

We also require that both the mass and radius of 
the exomoon can be considered independently detected. 
This means that the posteriors of both Ms/Mp and 
Rs / Rp must not be converged at zero in the planet-with- 
moon fits. These battery of tests form our requirements 
for an exomoon to be considered plausible and so far we 
have described four basic detection criteria, summarized 



Bl Improved evidence of planet-with-moon fits at > 
4tT confidence 

B2 Planet-with-moon evidences indicate both a mass 
and radius preference 

B3 Parameter posteriors are physical (e.g. pp, p^) 

B4 Mass and radius of moon converge away from zero 

Should these tests be passed or perhaps a candidate 
only marginally fails some of these criteria, we may con- 
sider further investigation. We discuss here three quick 
general follow-up criteria which can be implemented. 
Since this paper's survey only uses Ql-9 data, subse- 
quent Kepler data may be treated as follow-up photom- 
etry. One simple check then is that all four basic crite- 
ria are satisfied when the new data is included and the 
model refitted. Further, the significance of the moon 
candidate should be enhanced by the new data and yield 
broadly the same set of parameters. Without even fit- 
ting the new data, another simple check is to extrapolate 
the best-fit light curve model from the moon hypothesis 
into the times of the new observations and compare (in 
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a sense) the "predictive power" of the moon model 
relative to a simple planet-only model. We summarize 
these general follow-up detection criteria below: 

Fl All four basic criteria are still satisfied when new 
data is included 

F2 The predicitive power of the moon model is supe- 
rior (or at least equivalent) to that of a planet-only 
model 

F3 A consistent and statistically enhanced signal is re- 
covered with the inclusion of more data 

Even after passing all of these tests, the candidate 
should not be blindly accepted as a confirmed detection. 
Candidate specific tests and follow-up may be needed 
too, if for example the star shows rotational modulations 
and the candidate moon exhibits mutual events (which 
may in fact be star spot crossings). Exploration of per- 
turbing planet solutions causing TTVs/TDVs may be 
needed, as in the case of INesvornv et al.l (j2012D . Target 
specific tests will be investigated appropriately should 
the need arise. 

4.6. Excluded Moons 

In cases of null-detections, one of the aims of the HEK 
project is to provide limits on what moons can be ex- 
cluded. There are many possible choices for which pa- 
rameters we provide excluded limits. Two terms of par- 
ticular interest are the mass ratio (Mg/Mp) and the ra- 
dius ratio (Rs/Rp)- In general, we find the radius ra- 
tio to be untrustworthy due to the effects of starspots 
and time-correlated noise. Further, the parameter is 
positively-biased since we impose a likelihood penalty for 
high-density (i.e. low-radius) moon solutions. Therefore, 
we opt to provide upper limits on the mass ratio Ms/Mp 
for each null-detection. 

The excluded limits must be understood in terms of 
the adopted priors. For example, when we posit that 
there is a 95% probability of there existing no moon of 
Ms/Mp > X, the statement is only meaningful when 
combined with the adopted priors e.g. uniform prior in 
o-sp/Rp and Jeffrey's prior in P5. This subtle point is 
important when interpreting the 95% and 3 a quoted lim- 
its. However, we also make available the full posteriors 
revealing all relevant cross-correlations for those wishing 
to investigate the frequency of moons in more detail. 

One caveat with the provided upper limits is that 
MultiNest may have located a spurious signal and spu- 
rious detections cannot be used to define upper limits on 
Ms/Mp. Spurious detections occur because the code at- 
tempts to locate the best modes which explain the data 
i.e. the best model fits. In many cases, the solution 
can be dismissed using some of the detection criteria al- 
ready discussed but the derived posteriors of Ms/Mp 
and/or Rs/Rp still converge to non-zero values due to 
perhaps time-correlated noise or starspots. Simply tak- 
ing the 95% quantile of these posteriors does not techni- 
cally translate to an excluded upper limit estimate. In- 
deed these upper limits can approach unity if the model 
fit converged to a binary-planet solution, for example. In 
such spurious detection cases, all we can say for certain is 
that we are unable to detect an exomoon but we caution 



that meaningful upper limits on Ms/Mp and Rs/Rp is 
not guaranteed for each system analyzed. 

In addition to providing limits on Ms/Mp, we also 
compute limits for two more observable-centric defined 
terms, the moon-induced TTV and TDV amplitudes. 
Specifically, we calculate the quantiles of the distribu- 
tion of the root-mean-square (r.m.s.) amplitudes of th e 
two effects using the expressions fromlKiEEinJ (l2011bD . 
(Sttv and i5tdv- The TDV r.m.s. amplitude is defined as 
the sum of the TDV-V an d TDV- TIP effects, for reasons 
discussed in Chapter 6 of IKippin g (2011b). It should be 
stressed that excluding a certain moon-induced TTV am- 
plitude does not equate to excluding a TTV amplitude 
induced by other effects too (and similarly for TDV). 

5. RESULTS 
5.1. KOI-722.01 

5.1.1. Data selection 

After detrending with CoFiAM, the PA and PDC-MAP 
data were found to have a 1.7 a and 1.9 cr confidence of 
autocorrelation on a 30 minute timescale respectively and 
therefore both were acceptable (< 3cr). In general, we 
always prefer to use the raw data and so we opted for 
the PA data in all subsequent analysis of this system. 
Short-cadence data is available for quarter 9 and this 
data displaced the corresponding long-cadence quarter 
in our analysis. 

5.1.2. Planet-only fits 

When queried from MAST, the KIC effective temper- 
ature and surface gravity were reporte d as Tcb = 6133 K 
and log 5 = 4.628 (jBrown et al.ll2011[ ). Using these val- 
ues, we estimated quadratic limb darkening coefficients 
ui = 0.3694 and (ui -f- U2) — 0.6564 (as described in 
i j4.2[) . The initial two models we regressed were Vp and 
Vp,LD where the former uses the aforementioned limb 
darkening coefficients as fixed values and the latter al- 
lows the two coefficients to be free parameters. We find 
that logZ(Vp,LD) - logZ(Vp) = 0.05 ± 0.26 indicating 
essentially no preference between the two models. Given 
that the data is equally well-described by either theoret- 
ical or fitted coefficients, we opt for the theoretical limb 
darkening coefficients since they are more physically mo- 
tivated. 

KOI-722.01 has a period of Pp = 46.40630 ± 
0.00022 days (as determined by model Vp,ld) and ex- 
hibits 14 complete transits in Q1-Q9 from epochs -8 to 
-f8 (epochs -1-3, +4 and +7 are missing). As is typi- 
cal for all cases, logZ{J^p) > logZ(Vp) indicating that 
allowing for 14 independent baseline parameters is un- 
necessary relative to a single baseline term (thanks to 
the CoFiAM final stage normalizations). 

We find no evidence for TTVs in KOI-722.01, with 
logZ(J"TTv) - logZ(J'p) = -27.82 ± 0.18, which is for- 
mally an 7.2(7 preference for a static model over a TTV 
model. The timing precision on the 14 transits ranged 
fro m 5- 17 minutes (see Tabled]). The TTVs, shown in 
Fig. 2(a) show no clear pattern and exhibit a standard 



deviation of i5ttv = 12.1 minutes and Xttv ~ 21.6 for 
14-2 degrees of freedom. 

The TTV-hTDV model fit, Vy, finds consistent transit 
times with those derived by model J^ttv- We also find no 
clear pattern or excessive scatter in the TDVs, visible in 
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Fig. 2(a) The standard deviation of the TDVs is found 
to be (5xDV — 20.1 minutes and we determine Xtdv ~ 
15.7 for 14-1 degrees of freedom. 

5.1.3. Moon fits 

A planet-with-moon fit, J-s, is preferable to a planet- 
only fit at a formally high significance level (7.5 cr), pass- 
ing detection criteriopn Bl (see Table EJ. KOI-722.01 
fails detection criterion B2 though, since a zero-mass 
moon model fit yields a higher Bayesian evidence at 2.8 a 
confidence. 

Investigating further, one finds the parameter posteri- 
ors to be ostensibly unphysical. Physical parameters of 
the candidate solution may be estimated by combining 
our posteriors with the estimated stellar parameters of 
B12 (M^ = 1.08 Mq and = OMRq). The moon is 
found to lie at a highly inclined orbit and exhibit phys- 
ically consistent parameters of Ms = LlSlgygM© and 
Rs = 0.9061o;q55 i?©. In contrast, the planet has a 



very high density with parameters Mp = 245_y5 M© and 

Rp — 2.03llQ 04y i?©. We consider these parameters to 
be likely unphysical and thus KOT722.01 fails detection 
criterion B3. Finally, the posterior for Ms/Mp does not 
converge away from zero and thus no mass signal can 
be considered to be detected, failing criterion B4. This 
is also consistent with the fact that a zero-mass moon 
model gave a higher Bayesian evidence. We therefore 
conclude the model fit of J"s is an exomoon false-positive 
and no convincing evidence for a satellite around KOT 
722.01 exists in Ql-9. 

The origin of the false-positive is unclear as the max- 
imum a-posteriori fit from J-g reveals auxiliary transits 
(see Fig. which cannot be caused by starspot cross- 
ings. In this case, we consider time-correlated noise to 
be the most likely explanation. The Ms/Mp posterior 
converges on zero and we find the 95% quantile to be 
Ms/Mp < 0. 016 to the 3cr quantile is Ms/Mp < 0.031 
(see Fig. |4(a)| ). Our final system parameters are provided 
in Table IT"^ 

5.1.4. Summary 

We find no compelling evidence for an exomoon around 
KOT722.01 and estimate that Ms/Mp < 0.016 to 95% 
confidence. This assessment is based on the fact the sys- 
tem fails the basic detection criteria B2, B3 and B4 (see 

5.2. KOI-365.01 

5.2.1. Data selection 

After detrending with CoFiAM, the PA and PDC-MAP 
data were found to have a 2.2 cr and 3.1 tr confidence 
of autocorrelation on a 30 minute timescale respectively. 
Since only the PA data is considered acceptable (< 3cr), 
it will be utilized throughout the analysis that follows. 
Short-cadence data is available for quarters 3, 7, 8 and 9 
and this data displaced the corresponding long-cadence 
quarters in our analysis. 

5.2.2. Planet-only fits 

When queried from MAST, the KIC effective temper- 
ature and surface gravity were reported as Teff = 5389 K 



and log g = 4.570 (iBrown et al.ll2011[ ). Using these val- 
ues, we estimated quadratic limb darkening coefficients 
ui = 0.5016 and {U1+U2) = 0.7064. The initial two mod- 
els we regressed were Vp and Vp.ld, where the former 
uses the theoretical limb darkening coefficients as fixed 
values and the latter allows the two coefficients to be free 
parameters. We find that logZ(VpxD) — log-Z(Vp) = 
+5.34 ± 0.24 indicating that our limb darkening coeffi- 
cients may not be optimal. The maximum a-posteriori 
limb darkening coefficients from the Vp.ld model fit were 
ui = 0.5234 and (ui -f U2) = 0.6145, which we chose to 
treat as fixed terms in all subsequent model fits. 

KOI-365.01 has a period of Pp = (81.73766 ± 
0.00014) days (as determined by model Vp,ld) and ex- 
hibits 9 transits from Q1-Q9. As is typical for all cases, 
logZ(J^p) > logZ(Vp) indicating that allowing for 9 in- 
dependent baseline parameters is unnecessary relative to 
a single baseline term. 

We find no evidence for TTVs in KOT365.01, with 
logZ(J"TTv) - logZ(J"p) = -36.98 ± 0.17, which is for- 
mally an 8.3 a preference for a static model over a TTV 
model. The timing precision on the 9 transits ranged 
from 1.3 to 2. 8 minu tes and yields a fiat TTV profile, as 
shown in Fig. 2(b) We calculate a standard deviation 
of (5ttv ~ 1-5 minutes and x^tv ~ 4. 2 for 9-2 degrees of 
freedom. 

The TTV-FTDV model fit, Vy, finds consistent transit 
times with those derived by model J-ttv- No clear pat- 
tern orexcessive scatter is visible in the data, shown in 
Fig. 2(b) We therefore conclude there is no evidence for 
TTVsotTDVs for KOI-365.01. The standard deviation 
of the TDVs is found to be 5tdv — 3.3 minutes and we 
determine Xtdv ~ ^-^ degrees of freedom. 

5.2.3. Moon fits 

A planet-with-moon fit, J-s, is disfavored relative to 
a planet-only fit at A(logZ) = -2.7 ± 0.2 (see Tabled 
and so detection criterion Bl is not satisfied. The system 
also fails detection criterion B2 since a zero-mass moon 
yields an improved Bayesian evidence relative to a moon 
with finite mass. 

We may combine the posteriors from J^g with the stel- 
lar parameters derived by B12 (A/^ = 0.99 M© and 
i?* = 0.86 i?©) to obtain physical parameters for the 
planet-moon candidate system. We find that the planet 
has an unusually low density with Mp — 0.53j^Q'33 



and Rp — 2.1Y^^^^T^Rt^. The moon also has a low den- 



sity with parameters Ms — O.O6IQ O4M0 and Rs — 

0.65tQ 2i i?©- Whilst these parameters are somewhat 
extreme they are not implausible, but at best detec- 
tion criterion B3 can be considered unclear. Finally, the 
Ms/Mp posterior fails to converge away from zero, fail- 
ing detection criterion B4, and consistent with the fact 
that model J-"s,mo is favored over model J"s- We there- 
fore conclude the model fits of J-s and J^s,eP represent an 
exomoon false-positive and no convincing evidence for a 
satellite around KOI-365.01 exists in Ql-9. 

The maximum a-postcriori model fit of J^g (shown in 
Fig [SJ does not exhibit any clear auxiliary or mutual 
events, despite Rs/Rp converging away from zero. The 
reason for this becomes clear when one notes that the 
semi-major axis of the moon's orbit around the planet 
converges to (asp/Rp) = 3.7~!io'^, which shows that the 
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Fig. 2. — TTVs (black circles) and TDVs (gray squares) from model fit Vv of the seven KOls investigated here. TTVs computed relative 
to the maximum a-posteriori ephemeris derived by model J-p. TDVs computed relative to the median duration from model J-p. After the 
title of each sub-figure, we provide the confidence (in sigmas) of a TTV model over a static model in parentheses for each KOI (negative 
values imply a static model is favored). 



14 



Kipping et al. 



TABLE 4 

Transit times and durations for KOI-722.01. The model used to calculate the supplied values is provided in parentheses next to each 

column heading. BJDtj'^C times offset by 2, 400, 000 days. 



Epoch 


T [BJDuTc] (^TTV) 


TTV [mins] (J'ttv) 


T [BJDuTc] (Vv) 


TTV [mills] (Vv) 


T [mins] (Vv) 


TDV [mins] (Vy) 


-8 


54979.58041;;;;;;;^; 


-4.2 ±8.3 


54979.5798i;!;;!'^^» 


-12.0 ±7.0 


396I20 


-10.3 ±9.4 


-7 


55025.98891;;:;;;!^^ 


-2.2 ±9.9 


55025.992llg;'(!'(!^4 


-3.4 ±9.7 


4041^^ 


-6± 13 


-6 


55072.382+;;-;!i^ 


-22 ± 17 


55072.384+2'°^^ 


-24 ± 17 


466+^? 


+25 ± 26 


-5 


55ii8.8oioi;;;;S 


-4.9 ±8.0 


55ii8.8io+;!';!-:|^ 


+4 + 33 


464ljf 


+23 + 42 


-4 


55165.22111;;;";!;]^ 


+13.8 ± 5.4 


55i65.2i98i;!;;!;!« 


+9.0 ± 5.9 


406i;^ 


-5.8 + 7.4 


-3 


55211.61891;!;;!;!?;; 


+ 1 ± 11 


55211.63251;!;;!";^^ 


+18 + 15 


342+g 


-37 + 20 


-2 


55258.0267+0 ;!;!-g2 
55304.4253l'^-;»;!f, 


+1.7 ±8.4 


55258.0324i;!';!;!;!f, 


+9 + 14 


489+3f 


+36 + 21 


-1 


-10.4 ±8.3 


55304.42221!!;;;;;^? 


-15 + 10 


457li 


+20 + 13 


+0 


5535o.8409i;!;!;!« 


+2.0 + 6.9 


5535o.84i5i;!;;!;!;» 


3.5 + 8.2 


402118 


-8+ 11 


+1 


55397.246oi;;;;!;!'^^ 


-0.9 ±7.8 


55397.2466i;!';!o57 


1.7 + 8.8 


3981^;; 


-10 + 12 


+2 


55443.6628i;;;;!;!g 


+ 13.3 + 9.9 


55443.6668+;!-'2°?| 


+22 + 10 


39312S 


-12 + 14 


+5 


55582.8744l'(;;°°^? 


-0.3 + 8.4 


55582.874ll«;oo?^ 


+4.5 + 9.7 


4281^^ 


+5 + 13 


+6 


55629.28OOIH0S 


-2.4 ±5.9 


55629.2870lo gj74 


+14 ± 14 


448l« 


+15 + 17 


+8 


55722.0752l°;°°?5 


-30 ± 10 


55722.0784lo;°J24 


-17 ±15 


448lg 


+16 + 19 
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Fig. 3. — From left-to-right then top-to-bottom we show the chronological sequence of transits of KOI-722.01. The first 8 panels show 
the Ql-9 data and the maximum a-posteriori light curve fit of a planet-only model (gray line) and a moon model (red line). Note that the 
figure temporally zooms-in on the transits of the planet and the candidate moon. 
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Fig. 4. — Marginalized posterior distributions of Ms/Mp from model J-g for the seven KO Is investi gated here, except KOI-303.01 and 
KOI-1472.01 for which we use models J^s,RO and J^s,l3 respectively, for reasons discussed in i|5.6Rj5^ After the title of each sub-figure, 
we provide the 95% quantile on Ms/Mp in parentheses for each KOI. 
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TABLE 5 

Bayesian evidences of various fits for KOI-722.01. A description 
of the different models can be found in 



Model, M \ogZ{M) Mi - M2 

= logZiMi) - logZ(M2) 



Planet only fits. . . 
Vp 38458.54 ±0.18 

Vp,LD 38458.60 ±0.19 

Tp 38547.22 ±0.11 



Vp,LD - Vp = (+0.05 ± 0.26) 
- Vp = (+88.68 ± 0.21) 



Planet with timing variations fits... 

J^TTV 38519.40 + 0.14 J"ttv - -^P = (-27.82 ± 0.18) 

Vy 38283.44 ± 0.29 Vy - Vp = (-175.10 ± 0.34) 



Planet with moon fits... 
J^S 38577.91 ±0.12 

-^S,MO 38583.07 ±0.11 

•^330 38549.42 ±0.11 



J-g - J-p = (+30.68 ± 0.16) 
J^S.MO -^S = (+5.16 + 0.17) 
J^S.RO -^S = (-28.49 ±0.16) 



TABLE 6 

System parameters for KOI-722.01 from model Vp^LDj except for 
Mg/Mp which is derived from model Jfg. * indicates that a 
parameter was fixed. 



Parameter 



Value 



Derived parameters 
Pp [days] 
TO [BJDuTC 

Rp/R* 



46 40630+° °°°2^ 

2455350.8380lo 00I2 
0.02182«r?r 



b 

(a/R.) 
i [deg] 
Pi, [gcm-3 
T [hours] 

Ml 

(ui + U2) 



7+0.24 
-0.21 
1+4.4 
-9.0 



0.47 
46.1 

89.41+°;^? 




+0.27 
-0.41 
o lQrr+0.066 
o.±»o_o 058 
7+0.31 
-0.18 
^0.18 
-0.16 



0.273 

0.30^ 



Physical parameters. 
M* [i?0] 
R* [Rq\ 
Rp [R®] 
Ms/Mp 
Sttv [mins] 
(5tdv [mins] 



1.08* 

0. 83* 

1. »/o_0 039 



< 0.016 [95% confidence] 

< 4.0 (95% confidence) 

< 3.4 (95% confidence) 



moon is in very close proximity to the planet. In such 
a case, the planet and moon appear almost on-top of 
one another and thus virtually no light curve distortion 
is visible. This one way in which a fitting algorithm 
can essentially "hide a moon" in the fits and such fits 
are always suspicious. The close proximity of the moon 
leads t o an absence of TTVs too since TTVs scale as 
Mgas ('Kippingj l2009a|) and so the solution also masks 
the exomoon mass. 

Unlike TTVs, TDVs scale as Ms/y/os and so one 
might expect that the exomoon mass could not be hidden 
from TDVs too. However, the period of the moon solu- 
tion is also short at = 1.9ti;3 d ^^^^ ^^^^ causes a prob- 
lem for TDV inference. As noted in lKippin3 ()2011bD . 
traditional TDV theory breaks down if the moon accel- 
erates/decelerates significantly during a transit duration 
and thus the theory only holds for P5 ^ T14. In this 
case, Ps ~ 2 days and T14 '-^ 0.3 days meaning that that 
TDVs will not necessarily be present for even massive 
exomoons. Due to these points, close-binary moon so- 
lutions tend to yield less useful constraints on excluded 



mass and radius ratios. In this case, the moon fit yields 
a worse Bayesian evidence than a planet-fit and thus we 
know it is not the favored model, yet the upper lim- 
its are not as reliable as with KOT722.01 for example. 
The Ms I Mp posterior converges on zero and we find the 
95% quantile to be Ms/Mp < 0.69 and the Scr quan- 
tile is Ms/Mp < 0.97 (see Fig. |4 (b )|. Our final system 
parameters are provided in Table |9| 

5.2.4. Summary 

We find no compelling evidence for an exomoon around 
KOT365.01 and estimate that Ms/Mp < 0.69 to 95% 
confidence. This assessment is based on the fact the sys- 
tem fails the basic detection criteria Bl, B2 and B4 and 
B3 is considered marginal (see 

5.3. KOI- 174.01 

5.3.1. Data selection 

After detrending with CoFiAM, the PA and PDC-IMAP 
data were found to have a 1.5 tr and 0.6 ct confidence of 
autocorrelation on a 30 minute timescale respectively and 
therefore both were acceptable (< Scr). Since the PA 
data detrending shows no strong evidence of autocorre- 
lation, we opted to use this less-manipulated data in what 
follows. Short-cadence data is available for quarters 3, 4, 
5 and 6 and this data displaced the corresponding long- 
cadence quarters in our analysis. 

5.3.2. Planet-only fits 

When queried from IVIAST, the KIC effective temper- 
ature and surface g ravity were reporte d as Toff = 4654 K 
and log 3 = 4.538 JB rown et al.ll201l[ ). Using these val- 
ues, we estimated quadratic limb darkening coefficients 
ui = 0.6531 and (U1+U2) = 0.7415. The initial two mod- 
els we regressed were Vp and VpxD, where the former 
uses the theoretical limb darkening coefficients as fixed 
values and the latter allows the two coefficients to be free 
parameters. We find that logZ(Vp,LD) — log-Z(Vp) = 
+8.22 ± 0.22 suggesting that our limb darkening coeffi- 
cients could be improved. Strangely though, when we 
fix the limb darkening coefficients to the maximum a- 
posteriori values from the Vp^ld model fit, we actually 
obtain a worse Bayesian evidence than Vp. In light of 
this, we decided to continue with the theoretical coeffi- 
cients for this particular system. 

KOI-174.01 has a period of Pp = (56.35439 ± 
0.00019) days (as determined by model Vpxd) and ex- 
hibits full 9 transits from Q1-Q9. As is typical for all 
cases, logZ(J"p) > logZ(Vp) indicating that allowing 
for 9 independent baseline parameters is unnecessary rel- 
ative to a single baseline term. 

We find no evidence for TTVs in KOT174.01, with 
logZ(J"TTv) - logZ(J'p) = -22.65 ± 0.16, which is for- 
mally an 6.4 CT preference for a static model over a TTV 
model. The timing precision on the 9 transits ranged 
from 2.9 to 5. 9 min utes and yields a flat TTV profile, as 
We calculate a standard deviation of 
^2 - 13.0 for 9-2 degrees of 



shown in Fig. 2(c) 



(5ttv = 5.1 minutes and Xttv 
freedom. 

The TTV+TDV model fit, Vy, finds consistent transit 
times with those derived by model J-"ttv- The data show 



no clear pattern or excessive scatter, visible in Fig. 2(c) 
We therefore conclude there is no evidence for TTVs or 
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TABLE 7 

Transit times and durations for KOI-365.01. The model used to calculate the supplied values is provided in parentheses next to each 

column heading. BJDuTC times offset by 2, 400, 000 days. 



Epoch r [BJDuTc] (^TTv) TTV [mins] (J"ttv) r [BJDutc] (Vv) TTV [mins] (Vy) T [mins] (Vv) TDV [mins] (Vy) 



-5 


54962. 9413+00021 


+1.3 ± 2.9 


54962.9446+i;'^o^^ 


+6.1 ± 7.2 


-4 


55044.6779t«:ooi^ 


-0.2 ± 2.1 


55044.6775j:0-:ooil 


-0.6 ±2.1 


-3 


55126.4156+0 0011 


+0.0 ± 1.7 


55126.4160l«;»il 


+0.7 + 2.4 


-2 


55208.15191^0^1 


-1.9 ± 2.0 


55208.15211000^5 


-1.4 ±2.1 


-1 


55289. 89261o:oti5 


+2.4 ±2.3 


55289.8926lo ools 


+2.8 ± 2.1 


+1 


55453. 3657l|j:ggJ^ 


-0.7 ±2.7 


55453.3660to ooig 


+0.3 + 2.3 


+2 


55535. 1034+[5-{3^i2 


-0.5 ± 1.7 


55535. 1038^Q qq]^2 


+0.5 + 2.0 


+3 


55616.84083i«;r^^ 


-0.8 ± 1.4 


55616.84088lo;ooo93 


-0.1 + 1.4 


+4 


55698. 5805l:o:ooi2 


+2.1 ± 1.7 


55698. 5806lO;Ooi3 


+3.0 + 1.9 



402: 
382.5! 
386.4j 
391.1 



FTS — 

-14 
r+S.O 
^-5.3 
^+i.'o 
-6.3 
„+4.7 
3-5.0 

398.7l*-9 
397.ltt:^ 
387.3l*:0 
390.6l2g 

39i.oit;5 



+5.8 + 8.3 
-3.8 + 2.6 
-1.8 + 2.7 
+0.8 + 2.4 
+4.3 + 2.4 
+3.5 + 2.7 
-1.4 + 2.2 
+0.3 + 1.5 
+0.5 + 2.1 









/ * ' 

* * 



54-962.6 54 962.B 54963.0 54963.2 54963.4 




55044.2 55044.4 55044.6 55044.B 55045.0 



55126.0 55126.2 55126.4 55126.6 55126.B 





55207. B 5520B.O 5520B.2 5520B.4 55208.6 



552B9.4 552B9.6 55269.8 55290.0 55290.2 



55453.0 55453.2 55453.4 55453.6 55453.B 




55534.B 55535.0 55535.2 55535.4 55535.6 





55G1G.4 55616.E 55G1E.B 55617.0 



55E9a.2 556994 55699.6 55699.9 



Fig. 5. — From left-to-right then top-to-bottom we show the chronological sequence of transits of KOI-365.01. The first 8 panels show 
the Ql-9 data and the maximum a-posteriori light curve fit of a planet-only model (gray line) and a moon model (red line). No clear 
distortions due to the moon are visible in this close-binary type solution. 



TDVs for KOI-174.01. The standard deviation of the 
TDVs is found to be (5tdv — 6.7 minutes and we deter- 
mine XxDV ~ ^^-"^ ^'^^ degrees of freedom. 

5.3.3. Moon fits 

A planet-with-moon fit, J-s, is slightly favored rela- 
tive to a planet-only fit at A(logZ) — +1.5 ± 0.1, but 
does not satisfy detection criterion Bl. The system also 
fails detection criterion B2 since a zero-mass moon yields 
an improved Bayesian evidence relative to a moon with 
finite mass. 



Inspection of the posteriors from model J-s and us- 
ing the stellar parameters of B12 (M^, = 0.80 Mq and 
i?* ~ 0.80 M0) reveals a set of broadly unphysical 
parameters. Most notably, the planet has an unusu- 
ally low density with Mp ^ 0.27to^lMg) and Rp = 
2.28lg'27 i?e- "^^^ satellite also has some odd parame- 
ters with Ms = 0-08tom and Rs = 1.08t^;4o -^e- In 
general, we find this combination of masses and radii im- 
probable and consider that detection criterion B3 is not 
satisfied. Finally, the mass ratio Ms/Mp does not con- 
verge away from zero meaning detection criterion B4 is 
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TABLE 8 

Bayesian evidences of various fits for KOI-365.01. A description 
of the different models can be found in 



Model, M 


logZ{M) 


Ml - M2 
= \ogZ{Mx)- 




Planet only fits. . . 






Vp 


141526.62 ± 0.17 






Vp,LD 


141531.95 ± 0.17 


Vp.LD — Vp_ = 


(+5.34 ± 0.24) 


Vp,MAP 


141535.04 ± 0.17 


Vp,MAP — Vp 


= (+8.43 ± 0.24) 


^P 


141606.76 ± 0.10 


J"p — Vp,MAP 


= (+71.72 ±0.20) 


Planet with timing variations fits... 






141569.78 ± 0.14 


Tttv — J-'p = 


(-36.98 ± 0.17) 


Vv 


141394.37 ±0.26 


Vv — Vp^MAP 


= (-140.67 ±0.30) 


Planet with 


moon fits... 








141610.16 ± 0.12 


- -^P = (+ 


3.40 ±0.16) 




141611.69 ± 0.12 


■^S.MO — -^S = 


(+1.53 + 0.17) 




141602.99 ± 0.11 


-^S,RO - -^S = 


(-7.17 ±0.16) 



TABLE 9 

System parameters for KOI-365.01 from model Vp lDj except for 
Ms/Mp which is derived from model JFg. * indicates that a 
parameter was fixed. 



Parameter 



Value 



Derived parameters 
Pp [days] 
713 [BJDuTC 

Rp/R* 
b 

(a/R.) 
i [deg] 

T [hours] 

Ul 

{ui + U2) 



0.02364: 
0.47j 
84.5^ 



81 7Q7fifi+0 00014 
81.7376b_o 00014 

2455371. 62859t°E;gE;*^ 

+0.00088 
0.00058 

+0.18 
0.30 
+9.8 
-11.6 

89.68toig 
-I 7-1+0.66 

'^■"^-0.61 
R ^10+0.044 

Q4-+0.15 
652+° °33 



Physical parameters.. 



[Rq\ 

R* [R@] 
Rp [Rfs] 
Ms/Mp 
Sttv [mins] 
5tdv [mins] 



0.99* 
0.86* 

r, r,i 7 + 0.083 

' -0.054 

< 0.69 [95% confidence] 

< 3.3 (95% confidence) 

< 6.6 (95% confidence) 



also not satisfied. We therefore conclude that the model 
fit of J^g represents an exomoon false-positive and no con- 
vincing evidence for a satellite around KOT174.01 exists 
in Ql-9. 

The maximum a-posteriori model fit of J-'s (shown in 
Fig in]) does not exhibit any clear auxiliary or mutual 
events, despite Rs/Rp converging away from zero. The 
situation echoes that of KOI-365.01 and indeed both fits 
can be considered close-binary solutions. The planet- 
moon separation again converges to just a few planetary 
radii away {asp/Rp = 8.3^3;3), close enough that the 
moon appears essentially on-top of the planet in every 
transit. For the same reasons as described with KOT 
365.01, this close-binary solution results in poor con- 



straints on the exomoon mass. As shown in Fig. 4(c) 
the Ms/Mp posterior is unconverged yielding a 95% up 
per limit of Ms/Mp < 0.86 and a 3(7 upper limit of 
Ms/Mp < 0.99. Our final system parameters are pro- 
vided in Table [12] 



5.3.4. Summary 

We find no compelling evidence for an exomoon around 
KOT174.01 and estimate that Ms/Mp < 0.86 to 95% 
confidence. This assessment is based on the fact the sys- 
tem fails the basic detection criteria Bl, B2, B3 and B4 
(see S3]). 

5.4. KOI- 1472.01 

5.4.1. Data selection 

After detrending with CoFiAM, the PA and FDC-IVIAF 
data were both found to have a. 1.4 a confidence of auto- 
correlation on a 30 minute timescale and therefore both 
were acceptable (< 3 a). As before, we opt to use the FA 
data in what follows. No short-cadence data is available 
for this system and so long-cadence data only was used 
in what follows. 

5.4.2. Planet-only fits 

When queried from IVIAST, the KIC effective temper- 
ature and surface g ravity were reporte d as Teg = 5455 K 
and log (7 = 4.916 ([Brown et al.ll201l] ). Using these val- 
ues, we estimated quadratic limb darkening coefficients 
Ul = 0.4898 and (U1+U2) = 0.7037. The initial two mod- 
els we regressed were Vp and Vp,ld, where the former 
uses the theoretical limb darkening coefficients as fixed 
values and the latter allows the two coefficients to be free 
parameters. We find that logZ(Vp,LD) — log-Z^(Vp) = 
— 1.51 ± 0.22 suggesting that our limb darkening coeffi- 
cients are satisfactory. 

KOF1472.01 has a period of Pp = (85.35174 ± 
0.00020) days (as determined by model Vpxd) and ex- 
hibits 9 transits from Q1-Q9, with no transits absent 
from epoch -4 to +4. As is typical for all cases, 
logZ(J^p) > logZ(Vp) indicating that allowing for 9 in- 
dependent baseline parameters is unnecessary relative to 
a single baseline term. 

We find no evidence for TTVs in KOF1472.01, with 
logZ(J"TTv) - logZ(J-p) = -19.85 ± 0.17, which is for- 
mally an 6.0 a preference for a static model over a TTV 
model. The timing precision on the 9 transits ranged 
from 2.0 to 3.6 minutes (from model . Fttv ) and yields 



a flat TTV profile, as shown in Fig. |2(d) 
late a standard deviation of Jttv 

Xttv 

,2 



We calcu- 
4.6 minutes and 
Although the 



standard deviation of 
34.9 for 9-2 degrees of freedom. 
X" is somewhat excessive, the excess is not significant 
in light of the increased parameter volume required to 
produce these results i.e. the Bayesian evidence is lower 
than a static model. 

The TTV+TDV model fit, Vy, finds consistent transit 
times with those derived by model J-'ttv- As before, the 
dat a, show no clear pattern or excessive scatter, visible in 
Fig. 2(d) We therefore conclude there is presently no ev- 
idencelor TTVs or TDVs for KOF1472.01. The standard 
deviation of the TDVs is found to be Stuy = 6.0 minutes 
and we determine Xtdv = 34.2 for 9-1 degrees of free- 
dom. 

5.4.3. Moon fits 

A planet-with-moon fit, J\g, is favored relative to a 
planet-only fit at A(logZ) = +9.5 ± 0.2 (see Table [Ti]). 
The formal significance equates to 3.96 cr and can there- 
fore be considered to lie on the margin of satisfying de- 
tection criterion Bl. 



20 



Kipping et al. 



TABLE 10 

Transit times and durations for KOI-174-01. The model used to calculate the supplied values is provided in parentheses next to each 

column heading. BJDuTC times offset by 2, 400, 000 days. 



Epoch T [BJDuTc] (-^TTv) 



TTV [mins] (J"ttv) t [BJDutc] (Vy) TTV [mins] (Vv) T [mins] (Vv) TDV [mins] (Vy) 



-6 


54977.832ll'^;;jj;^^ 


-4.3 ±5.2 


54977.8295+^;^^* 


-8.0 ±6.5 


207+^7 


+8.5 + 8.5 


-5 


55034. 1934l0;°028 


+5.8 ±4.0 


55034. 1933lo 0029 


+5.7 ±4.2 


190.31^03^ 


+0.0 + 5.0 


-4 


55090. 5427+°oE52g 


-1.1 ± 3.8 


55090.5397lo 0030 


-5.5 ±4.7 


165+11 


-12.6 + 6.0 


-3 


55146.8945lO-;g;520 


-4.6 ±2.9 


55146. SdAStomll 


-4.3 ±4.9 


187.9t«-| 


-1.2 ±4.6 


-2 


55203.2535 q qq^q 


+2.2 ±5.9 


55203.2560lEj;{j|]g 


+5.8 ±5.5 


203^15 


+6.3 + 6.6 


+0 


55315. 9626log°24 


+3.1 ± 3.5 


55315.962lt00025 


+2.4 ±3.4 


198.7t^;g 


+4.2 + 3.9 


+2 


55428. 6705togo2l 


+2.2 ±2.7 


55428.6705looo2i 


+2.3 + 3.0 


I93.6tt:l 


+1.6 + 3.3 


+5 


55597. 7259tO-;«EJ|° 


-8.3 ±4.3 


55597.7263t«;°o^^ 


-7.8 ±4.8 


n9tll 


-5.9 ±6.0 


+6 


55654.0904tO;S 


+6.5 ±4.4 


55654.090i1E;;°o^8 


+6.0 + 4.0 


^7S.6tl%° 


-5.9 ±4.9 





54977.4 54977.6 54977.B 5497B.0 5497B.2 



55D33.B 55034.0 55034.2 55034.4 55034.6 



55090.2 55090.4 55090.E 55a90.B 55091.0 




.4 55146.6 55146. B 55147.0 55147.2 55147.4 




i5202.B 55203.0 55203.2 55203.4 55203.6 



55315.6 55315. B 55316.0 55316.2 55316.4 









ii 





Fig. 6. — From left-to-right then top-to-bottom we .show the chronological sequence of transits of KOI-174.01. The first 8 panels show 
the Ql-9 data and the maximum a-posteriori light curve fit of a planet-only model (gray line) and a moon model (red line). No clear 
distortions due to the moon are visible in this close-binary type solution. 



The posteriors reveal a solution hitting against the 
lower boundary condition on the planet-moon separa- 
tion, yielding asp/Rp — ?>.2+]^^ ■ The proximity be- 
tween the two objects is particularly extreme when one 
notes that the fit also finds a substantial satellite size of 



Rs/Rp = 0.40: 



-0.16 



The solution is therefore consistent 



with a close-binary, as we also found for KOI-365.01 and 
KOI-174.01. Due to the low density of the planet solution 
relative to that of the satellite, the posteriors indicate a 
Roche limit well-inside the planet and thus we cannot 
exclude the satellite with a tidal disruption argument. 



The stellar parameters for the host star are estimated 
in B12 as M^, = 0.93 and i?* = 0.56 i?0, from 
which we can estimate physical parameters for the can- 
didate planet-moon system. From model J"sj we deter- 
mine planetary parameters of Mp = IS.G^^^'g M0 and 
Rp = 3.501q'2i R® indicating a low density of pp = 
O.STlJ'gg gcm~'^, making the planet a small Neptune- 
type planet. The satellite returns Ms = O.SI^^q^ and 



Rs — 1.40to'46 suggesting a Super-Earth type moon. 



The satellite's period is Ps = 0.89t^;4^ days and thus 
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TABLE 11 

Bayesian evidences of various fits for KOI-174-01- A description 
of the different models can be found in ^4-l\ 



Model, M 


logZ{M) 


Ml - M2 
= logZ{Mi)- 


- logZ{M2) 


Planet only fits. . . 






Vp 


68914.54 ±0.16 






Vp,LD 


68922.76 ±0.16 


Vp.LD — Vp_ = 


(+8.22 ±0.22) 


Vp,MAP 


68914.48 ±0.16 


Vp,MAP — Vp 


= (-0.06 ±0.22) 


^P 


68987.21 ±0.10 


J'p — Vp,MAP 


= (+72.67 ±0.19) 


Planet with timing variations fits... 






68964.57 ±0.13 


Jttv -^P = 


(-22.65 ±0.16) 


Vv 


68825.27 ±0.21 


Vv — Vp^MAP 


= (-89.27 + 0.26) 


Planet with 


moon fits... 






J's 


68988.75 ±0.10 


J^S - J-p = (+1.54 ±0.14) 




68991.09 ±0.10 


-^S,MO ~ -^S = 


(+2.34 + 0.14) 




68979.64 ±0.11 




(-9.12 + 0.15) 



TABLE 12 

System parameters for KOI-174.01 from model Vp lDj except for 
Mg/Mp which is derived from model Tq. * indicates that a 
parameter was fixed. 



Parameter 



Value 



Derived parameters.. 
Pp [days] 
To [BJDuTC 

Rp/R* 
b 

(a/R^) 
i [deg] 
p* [gcm-^] 
T [hours] 

Ul 

{ui + U2) 



56 3543Q+0 00018 

2455315.96099^0 °g°75 

0.02908to°r? ■ 
0.46l«:- 
119.5t^-6 

89.78lin82 
10.2^2 3 

q iQcr+b.066 
■J-r»o_o 058 

u.au_Q 26 



Physical parameters.. 



AL [Rq] 
R* [Re] 
Rp [-Re] 
Ms/Mp 
5ttv [mins] 
5tdv [mins] 



0.80* 
0.80* 

O t:c!7+0.083 

< 0.86 [95% confidence] 

< 6.5 (95% confidence) 

< 5.8 (95% confidence) 



low enough (relative to the transit duration of 0.26 days) 
to conveniently mask TTV and TDV effects. As with 
the previous two close-binary solutions, radius-effects are 
also difficult to identify in the maximum a-posteriori fit 
shown in Fig. [71 Although the satellite candidate is just 
a few planetary radii away from its host, the posteriors 
are broadly physical and we therefore conclude that de- 
tection criterion B3 is satisfied. 

Out of all of the moon fits attempted, the highest 
Bayesian evidence comes from J^s,RO zero-radius moon 
model) as shown in Table [TH meaning that detection cri- 
terion B2 is failed. We note that the Ms/Mp ratio is not 
well- converged in the J^s model, but does appear diver- 
gent from zero, broadly satisfying criterion B4. 

In conclusion, KOI-1472.01 passes Bl, B3 and B4 but 
fails criterion B2. We therefore consider that more data 
may resolve whether the moon hypothesis is plausible or 
not. 

5.4.4. Predictive power of the moon model 



From a perspective, the maximum a-posteriori re- 
alization from model J-"s is naturally lower than that of 
J^P where we find 2181.01 versus 2248.48 respectively, for 
1495 data points spanning Ql-9. At the time of writing, 
QlO-13 had recently become available and this data may 
be used as a test between the planet and planet-with- 
moon hypotheses. If our moon model is genuine, then 
extrapolating the model into QlO-13 should yield a bet- 
ter prediction (in a sense) than the simple planet-only 
model. We downloaded and detrended this data accord- 
ingly using CoFiAM and the PA time series, which covers 5 
new transits and 815 new LC data points. Given that we 
got a improvement of 67.5 over Ql-9 (1495 points), we 
might expect the moon model to yield an improvement 
of - 30 in QlO-13. 

We extrapolated the maximum a-posteriori realization 
of model J"s into QlO-13 and find = 1222.85 for 
815 data points. Repeating the process for yields 
= 1196.54. We therefore find that the moon model 
has substantially worse predictive power than a simple 
planet-only model. Since the maximum a-posteriori re- 
alization is just a single realization, we decided to ex- 
trapolate 50 realizations randomly drawn from the joint 
posteriors of J^s in order to account for the parameter 
uncertainties (shown in Fig. [T])- From these 50 realiza- 
tions, we compute 50 x^ values from which we take the 
mean to be 1233 ± 22. Repeating the process for J^p 
yields x^ = 1201.7 ± 8.5, again supporting the planet- 
only model over the planet-with-moon model. 

These results clearly show that the planet-only model 
has superior predictive power to that of the moon model, 
thus the moon hypothesis fails detection criterion F2. 

5.4.5. Refitting the updated data set 

The Ms/Mp posterior derived from Ql-9 is converged 
off-zero due to a spurious detection. In order to place 
constraints on this this parameter, we require a posterior 
derived from a null detection instead. To this end, we 
decided to re-fit the Ql-13 data with an updated moon 
model, J-s,i3- A secondary goal was to check whether 
detection criterion F3 could be satisfied by refitting i.e. 
the same signal is retrieved as before but with higher 
significance. 

The regression yielded a distinct solution with both the 
semi-major axis and period of the moon moving outwards 
by many sigma. Specifically, we find asp/Rp = 69^}^ 
and Ps = 29.91y^7^ days. The mass ratio posterior, 
Af 5 /Mp , also shows a distinct profile and now converges 
close to zero with Ms/Mp = 0.0117tQ;oo52, which is con- 
sistent with a null-detection. The fact the Ms/Mp pos- 
terior converges close to zero means that criterion Fl is 
not satisfied and the fact that the solution is distinct 
from that obtained by the Ql-9 data means that crite- 
rion F3 is also not satisfied. In summary then, all three 
follow-up criteria are not satisfied and on this basis we 
conclude that there is no evidence for an exomoon around 
KOI-1472.01. 

In Table [T5l we display our final system parameters for 
KOI-1472.01. Given that the Ql-9 fits yielded a spurious 
moon detection but the Ql-13 fits are consistent with a 
null-detection, we use the latter for our constraints on 
Ms /Mp. The posterior of Ms/Mp from J^s 13, shown in 
Fig. [4(d)| yields a 95% quantile of Ms/Mp < 0.037 and 
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Fig. 7. — From left-to-right then top-to-bottom we show the chronological sequence of transits of KOI-1472.01. The first 8 panels show 
the Ql-9 data and the maximum a-posteriori light curve fit of a planet-only model (gray line) and a moon model (red line). The last 5 
panels show the QlO-13 data with 50 extrapolations (50 different shadings used) of the moon model overlaid (parameters randomly drawn 
from the joint posteriors) , which exhibit poor predictive power. 



II. The Hunt for Exomoons 



23 



TABLE 13 

Transit times and durations for KOI-1472.01. The model used to calculate the supplied values is provided in parentheses next to each 

column heading. BJDuTC times offset by 2, 400, 000 days. 

Epoch T [BJDuTc] (.T^TTv) TTV [mins] (J'ttv) t [BJDutc] (Vy) TTV [mins] (Vv) T [mins] (Vv) TDV [mins] (Vy) 

^4 54994.0093+o oo{4 +5.1 ±2.0 54994.0093+i;;ooJ4 5.81 ± 2.0 38r9Tf^ 2.9 ± 2.3 

-3 55079.3573t°;°°26 -0.3 ± 3.6 55079.3530l!];""2o -6 ± 14 394^^^ 6 ± 16 

-2 55164.7086l!;;!',!',i| -0.9 ±2.1 55164.70861!^;'^^^^ -0.4 ±2.1 383.81^1 0.9 ± 2.5 

-1 55250.0537lo 001^ -10.5 ± 2.2 55250.0554lE;;5;5ji^ -7.5 ± 2.8 362.611;° -9.7 ±3.1 

+0 55335.4153lo;ooJ| +3.7 ± 2.2 55335. 4153lo;ool| 4.0 ± 2.2 386.8I5 J 2.4 ±2.5 

+1 55420.7624l{5;g^^| -3.1 + 2.3 55420.7626l!];ooi6 -2.5 ± 2.3 401.4l|;f 9.7 ±2.6 

+2 55506.1166l°;ogi4 +0.5 + 2.0 55506.1169l„;""i;5 1.2 + 2.2 373.111;^ -4.5 + 2.5 

+3 55591.4683l°;°!3ig +0.4 + 2.8 55591.4683l„;ooi9 0.5 + 2.9 387.7lg;4 2.8 + 3.2 

+4 55676.8215l[j;5^5^i| +2.6 + 2.1 55676.8216l55;'(]^J^ 2.8 + 2.0 373.51^;^ -4.3 + 2.3 
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TABLE 14 

Bayesian evidences of various fits for KOI- 1472. 01. A 
description of the different models can be found in ^4-l\ The 
"13" subscript denotes that Ql-13 data was used rather than 
Ql-9. 

Model, M \ogZ{M) Mi - M2 

=logZ{Mi)-logZ{M2) 

Planet only fits. . . 
Vp 9485.63 ±0.15 - 

Vp.LD 9484.12 ±0.15 

Tp 9541.89 ±0.10 



Vp,LD — 
JFp - Vp 



Vp = (-1.51 ±0.22) 
= (+56.26 ±0.18) 



Planet with timing variations fits... 

J^TTV 9522.04 ± 0.13 J^ttv " -^p = (-19.85 ± 0.17) 

Vy 9385.63 ± 0.21 Vy - Vp = (-100.00 ± 0.26) 



Planet with moon fits... 
J^S 9551.38 ±0.11 

-^S,MO 9550.19 ±0.11 

J"s'ro 9552.32 ±0.11 



J-g - jTp = (+9.49 ±0.15) 
^S,MO -^S = (-1.18 ±0.16) 
^S,RO -.^S = (+1.03 + 0.16) 



-^P,13 
-^S.13 



14782.31 ±0.10 
14810.85 ±0.13 



p. 13 



(+28.54 + 0.16) 



TABLE 15 

System parameters for KOI-1472.01 from model Vp^LD, except 
for Mg/Mp which is derived from model J^s,i3- * indicates that 
a parameter was fixed. 



Parameter 



Value 



Derived parameters... 
Pp [days] 
TO [BJDuTc] 

Rp/R* 
b 

(a/R.) 
i [deg] 
P* [gcm-3] 
t [hours] 
ui 

(ui + U2) 



Rc; Qtr-i 7/1+0.00020 

00.001 fi.QQQQJQ 

2455335.41282toog05^ 
06422+° '"'2* 
U.OZ_Q 33 
88.51^^7'^ 

89.66to?4 

l-OU_Q 53 

6 276+° °*! 
0.63t°-3i 
0.53 



-0.22 

+0.17 
-0.13 



Physical parameters. 



M. [Rq] 
R* [Rq] 

Ms/Mp 
<5ttv [mins] 
<5tdv [mins] 



0.93* 
0.56* 
3 

■^•"^-0.14 

< 0.037 (95% confidence) 

< 2.7 (95% confidence) 

< 2.3 (95% confidence) 



a 3ct quantile of Ms/Mp < 0.063. 

5.4.6. Summary 

We find no compelling evidence for an exomoon around 
KOI-1472.01 and estimate that Ms/Mp < 0.037 to 95% 
confidence. This assessment is based on the fact the sys- 
tem fails the basic detection criterion B2 as well as the 
follow-up criteria Fl, F2 and F3 (see ^4.51) . 

5.5. KOI- 1857.01 
5.5.1. Data selection 

After detrending with CoFiAM, the PA and PDC-MAP 
data were found to have a 2.3 cr and 2.3 tr confidence of 
autocorrelation on a 30 minute timescale respectively and 
therefore both were acceptable (< 3 cr). As with previous 
systems, we choose to use the PA data over the PDC- 
IMAP data as both are acceptable. No short-cadence data 



is available for this system and so long-cadence data only 
was used in what follows. 

5.5.2. Planet-only fits 

When queried from IVIAST, the KIC effective temper- 
ature and surface gravity were reporte d as Teff = 5619 K 
and logg = 4.527 (jBrown et al.l[20TTI ). Using these val- 
ues, we estimated quadratic limb darkening coefficients 
Ml = 0.4567 and (U1+U2) = 0.6910. The initial two mod- 
els we regressed were Vp and Vp^ld, where the former 
uses the theoretical limb darkening coefficients as fixed 
values and the latter allows the two coefficients to be free 
parameters. We find that logZ(Vp^LD) ~ log2^(Vp) = 
+0.58 ± 0.21 suggesting that our limb darkening co- 
efficients may not be optimal. We thus chose to set 
the coefficients to the maximum a-posteriori values of 
ui = 0.5409 and (ui + us) = 0.8872. 

KOI-1857.01 has a period of Pp = (88.64486 ± 
0.00075) days (as determined by model Vp,ld) and ex- 
hibits 8 transits from Q1-Q9 from epoch -4 to +4, ex- 
cept epoch -2 which is absent. As is typical for all cases, 
logZ(J^p) > log2^(Vp) indicating that allowing for 8 in- 
dependent baseline parameters is unnecessary relative to 
a single baseline term. 

We find no evidence for TTVs in KOI-1857.01, with 
logZ(J"TTv) - logZ(J"p) = -18.59 ± 0.15, which is for- 
mally an 5.8 tJ preference for a static model over a TTV 
model. The timing precision on the 8 transits ranged 
from 6.8 to 9.5 minutes (from model Tttv ) and yields 



a flat TTV profile, as shown in Fig. 2(e) We calcu- 
late a standard deviation of Stty = 11.0 minutes and 
Xttv = ^^-^ degrees of freedom. 

The TTV+TDV model fit, Vy, finds consistent tran- 
sit times with those derived by model J-ttv- We detect 
no clea r patt ern or excessive scatter in the data, visible 



in Fig. 2(e)[ We therefore conclude there is presently 
no evidence for TTVs or TDVs for KOI-1857.01. The 
standard deviation of the TDVs is found to be Stuv = 
18.6 minutes and we determine Xtdv ~ "^-"^ ^'^^ ^"^ 
grees of freedom. 

5.5.3. Moon fits 

A planet- with- moon fit, J-"s, is slightly favoured rel- 
ative to a planet-only fit at A(logZ) = +2.16 ± 0.14, 
or 1.57 tJ, meaning detection criterion Bl is not satisfied. 
The zero-mass moon model is also preferred meaning cri- 
terion B2 is not satisfied either. 

Despite failing Bl and B2, the fits yield broadly physi- 
cal parameters and the low confidence could be indicative 
of a low signal-to-noise moon embedded in the data. Us- 
ing = 0.94 Mq and = 0.85 i?© from B12, we de- 
termine Mp = 6.1 



Ms 



1+19.3 7.^ 
-4.6 

O.bAtl fgM^foT Rs 



for Rp = 2.0671°;°^^ i?e and 



0.: 



3+0.14 
'-0.16 



Rffi . The solution 



is also clearly not a close-binary with asp/Rp = 42^25 
and auxiliary transits driving the fit, as evident from 
Fig. [51 Finally, the mass and radius ratio posteriors show 
convergence away from zero meaning criteria B3 and B4 
are satisfied. 

5.5.4. Predictive power of the moon model 

At this stage, we considered KOI-1857.01 to be a po- 
tential candidate and further data may confirm/reject 
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TABLE 16 

Transit times and durations for KOI- 1857.01. The model used to calculate the supplied values is provided in parentheses next to each 

column heading, BJDnq^c times offset by 2, 400, 000 days. 



Epoch T [BJDuTc] (-^TTv) 
0+U.UU48 



TTV [mins] (J'ttv) r [BJDutc] (Vy) TTV [mins] (Vv) T 



mms 



(Vv) TDV [mins] (Vy) 



^+U.UUbU 
-0.0047 

55067. 1017tE;;°°^^ 
55244.4029j;o 
55333. 0430loo°gg 
55421.68961°;°°^^ 

55510.3446lo 0070 
55598.9899l°;°°g 

,+0.0088 



TTTT- 

46lli 
,+34 
-23 
n+23 
-22 
,+21 
-19 
h29 
-24 

4501?^ 

H28 
-53 



-4 

-3 

-1 
+0 
+1 
+2 
+3 
+4 



54978.4743 
55067.1001 
55244.4050 
55333.0409 
55421.6891 
55510.3429 
55598.9909 
55687.6201 



0.0047 
+0.0056 
-0.0056 
1+0.0047 
-0.0047 
,+0.0051 
■0.0051 
+0.0054 
-0.0050 
,+0.0061 
■0.0058 
1+0.0057 
-0.0059 
+0.0064 
0.0068 



+11.6 ± 6.9 
-16.2 ± 8.1 
+4.8 ±6.8 
-8.5 ±7.3 
-4.1 ± 7.5 
+8.4 ±8.6 
+12.5 ± 8.4 
-10.4 ±9.5 



54978.4744^ 



55687.6241" 



0.0080 



+7.1 ± 7.0 
-18± 11 
-0.2 ± 9.7 
-6.5 ± 8.8 
-3.6 ± 8.3 
+11.7 ±9.7 
+12.8 ± 7.6 
-2 ±12 



471"* 



493: 
450j 
47i:! 
507"^ 



528"' 



+8 ±9.7 
+3± 14 

+18 + 14 
-3 ± 11 
+7 + 10 

+25 + 13 
-3 + 10 

+36 ± 45 



'itt lull III ,i till ,ii 
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Fig. 8. — From left-to-right then top-to-bottom we show the chronological sequence of transits of KOI-1857.01. The first 8 panels show 
the Ql-9 data and the maximum a-posteriori light curve fit of a planet-only model (gray line) and a moon model (red line). The last 4 
panels show the QlO-13 data with 50 extrapolations (50 different shadings used) of the moon model overlaid (parameters randomly drawn 
from the joint posteriors), which do not exhibit significant predictive power. 
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the signal. We therefore detrended the PA data for KOI- 
1857 from QlO-13 covering 4 new transits, which had 
recently become available at the time of writing. Before 
we attempted to re-fit the updated data set, we extrap- 
olated the light curve model fit J^g into QlO-13 for both 
the maximum a-posteriori solution and 50 randomly sam- 
pled solutions from the joint posteriors, which are shown 
in Fig. |8l The maximum a-posteriori moon solution can 
be compared to the detrended data, where we compute 
= 744.8 for 621 data points. For comparison, we re- 
peated the process for the maximum a-postcriori planet- 
only model and compute — 743.6 i.e. nearly identical 
but slightly worse. 

This comparison suggests that the moon model has no 
significant predictive power, which is a major concern 
for accepting the moon hypothesis. We should expect a 
moon model to give a better by around ^10 based 
upon the improvement of the best fits on the Ql-9 data. 
For Ql-9, we found Ax^ — 30.3 between the two best 
fits for 1437 data points. The 50 realizations of light 
curve predictions, shown in Fig. |8l do not seem to show 
a convincing agreement with the data either. Indeed, the 
distribution of the values for these fits yields x^ = 
724 ± 12 whereas repeating the process for 50 predictions 
from Tp yields 719 ± 10, which again suggests that the 
planetary model actually has slightly better predictive 
power. 

5.5.5. Refitting the updated data set 

The moon model J^s exhibits approximately the same 
predictive power as a simple planet-only model (although 
technically slightly worse) . As the difference in predictive 
power is not substantial (unlike with what we had with 
KOI-1472.01), we considered that the moon hypothesis 
could not yet be excluded. To fully exploit the QlO-13 
data, we stitched the detrended data onto the previous 
Ql-9 time series and re-fitted the updated data with a 
planet-only model ( J-pis) and a planet-with-moon model 

(^S13)- 

We find that the significance of the moon solution has 
indeed been enhanced moving from 1.57 cr to 3.81 cr, al- 
though still below the threshold of detection criterion 
Bl. 

The new solution finds a much more massive planet 
than before with Mp = 82^11° for Rp = 
2.3llo.i4 -^e- The satellite appears broadly unchanged 
in terms of mass and radius with Ms = 0.59^^^'^ for 
Rs = 0.97lto;QgQi?©. Despite this, the orbital solution 

is distinct with asp/Rp = 78^42 i.e. the moon is at 
about twice the orbital separation than before. One sig- 
nificant difference is that with just Ql-9 data, both the 
planet-only fit and the planet-with-moon fit obtained a 
Pi, value which appeared consistent with the KIC spec- 
tral classification (late G). Specifically, we found — 
1-4^0 4 gcm"'^ from the moon fit and from 1.3lQ ygcm~'^ 
from the planet fit. However, the J^s.ia model pushes 
down to 0.531^:1^^ gem meaning that the star is now 
more consistent with an F-dwarf. The KIC catalogue is 
known to be inaccurate for K/M dwarfs but G/F dwarfs 
are quite reliably identified and thus this discrepancy is 
suspicious. It is possible the planet is eccentric but given 
the wide orbit of the moon solution, even a moderate ec- 



TABLE 17 

Bayesian evidences of various fits for KOI- 1857. 01. A 
description of the different models can be found in ^4jJ\ The 
"13" subscript denotes that Ql-13 data was used rather than 
Ql-9. 



Model, M 


logZ{M) 


Ml - M2 








= \ozZ{Mi)- 


- \ogZ{M2) 


Planet only fits... 






Vp 


10535.69 ±0.15 






VpxD 


10536.27 ±0.15 


Yp,ld — Vp_ = 


(+0.58 ±0.21) 


Vp,MAP 


10537.26 ±0.15 


Vp,MAP — Vp 


= (+1.57 + 0.21) 




10592.10 ±0.09 


_Fp _ Vp,MAP 


= (+54.84 + 0.18) 


Planet with timing variations fits... 






10573.51 ±0.12 


.T^TTV — .^P = 


(-18.59 + 0.15) 


Vv 


10445.32 ± 0.22 


Vv — Vp^MAP 


= (-91.94 + 0.27) 


Planet with 


moon fits... 








10594.26 ±0.10 


- ^P = (+ 


2.16 + 0.14) 




10602.03 ±0.11 


-^S,MO - -^S = 


(+7.77 + 0.15) 


■^S,RO 


10592.41 ±0.10 


-^S,RO - -^S = 


(-1.85 + 0.15) 


.^P,13 


15158.41 ±0.10 








15167.30 ±0.10 


■^S,13 — ■^P,13 


= (+8.89 + 0.15) 


.?^S,M0,13 


15177.36 ±0.12 


^SMOXi - -^3,13 = (+10.06 ± 0.16) 


.?^S,R0,13 


15156.57 ±0.10 


.?^S,R0,13 — .^S 


13 = (-10.73 + 0.15) 



centricity would mean that the moon was unstable. We 
consider that detection criterion B3 is marginally satis- 
fied. However, in these new fits the Ms/Mp posterior 
now appears well- converged on zero, failing criterion B4. 

With Ql-9, we found that the zero-mass moon model 
was preferred over the physical moon model with 
AlogZ — 7.77±0.15, thus failing detection criterion B2. 
This result is further supported by the Ql-13 analysis 
where the same comparison yields A log Z = 10.09±0.16. 
We therefore find that follow-up detection criteria Fl, 
F2 and F3 are all not satisfied and the candidate can 
be dismissed. There is therefore no evidence for an exo- 
moon around KOI-1857.01 from Ql-13. From the Ql-9 
analysis only, the Ms/Mp posterior converges off zero 
and indicates Ms/Mp < 0.49 to 95% confidence and 
Ms/Mp < 0.95 to 3 cr confidence. Including the QlO-13 
data eliminates the mass convergent solution and leads 
to much tighter constraints of Ms/Mp < 0.028 to 95% 
con fidenc e and Ms/Mp < 0.052 to 3cr confidence (see 
Fig.[4(i)). 



5.5.6. Summary 

We find no compelling evidence for an exomoon around 
KOI-1857.01 and estimate that Ms/Mp < 0.028 to 95% 
confidence. This assessment is based on the fact the sys- 
tem fails the basic detection criteria Bl and B2 as well 
as the follow-up criteria Fl, F2 and F3 (see ^4.5p . 

5.6. KOI-303.01 
5.6.1. Data selection 

After detrending with CoFiAM, the PA and PDC-MAP 
data were found to have a 2.0(7 and 2.6 cr confidence of 
autocorrelation on a 30 minute timescale respectively and 
therefore both were acceptable (< 3cr). The less ma- 
nipulated PA data is selected over the PDC-MAP data 
throughout. Short-cadence data is available for this sys- 
tem for quarters 6, 7 and 8 and this data displaced the 
corresponding long-cadence time series in what follows. 

5.6.2. Improved stellar parameters 
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TABLE 18 

System parameters for KOI-1857.01 from model Vp^LD, except 
for Mg/Mp which is derived from model J^s,l3- * indicates that 
a parameter was fixed. 



Parameter 



Value 



0.43^ 



Derived parameters.. 
Pp [days] 
TO [BJDUTC] 

Rp/R* 
b 

{a/R,) 
i [deg] 
P* [gcm~^] 
f [hours] 
ui 

(mi + U2) 



2455333.0474l°°°2o 

,+0.33 
-0.30 

89.70to 43 

7+0.40 
-0.77 
-+0.18 
-0.17 
,+0.32 
-0.32 
,+0.13 



1.27: 

7.55^ 
0.72^ 
0.81^ 



0.18 



Physical parameters. . . 

M* [Rq] 0.94* 

R^, [Rq] 0.85* 

Rp [i?e] 2.210t0;g|? 

Ms/Mp < 0.028 (95% confidence) 

<5ttv [mins] < 5.9 (95% confidence) 

<5tdv [mins] < 7.7 (95% confidence) 



During the analysis of this candidate, we were able to 
acquire two high resolution spectra of KOI-303.01 using 
the ARCES spectrograph on the 3.5 m Astrophysical Re- 
search Consortium Telescope at the Apache Point Obser- 
vatory on 31 August 2012. We used a 1'.'6 x 3'.'2 slit with 
an exposure time of 20 minutes yielding a signal-to-noise 
of 25 per pixel with a resolution of 31,500. The echelle 
spectra were extracted using standard IRAF tools. 

We us ed the Stellar Param eter Classification (SPC) 
method (jBuchhave et al.l 120121 ) to derive the stellar at- 
mosphere parameters. SPC cross-correlates the observed 
spectrum against a grid of synthetic spectra drawn from 
a library calcu lated by John Laird using Kurucz models 
(|Kurucdll992[ ) . The synthetic spectra cover a window of 
300 A centered near the gravity-sensitive IVIgb features 
and have a spacing of 250 K in effective temperature, 
0.5 dex in gravity, 0.5 dex in metallicity and lkms~^ 
in rotational velocity. To derive the precise stellar pa- 
rameters between the grid points, the normalized cross- 
correlation peaks were fitted with a three dimensional 
polynomial as a function of effective temperature, sur- 
face gravity and metallicity. This procedure was carried 
out for different rotational velocities and the final stellar 
parameters were determined by a weighted mean of the 
values from the spectral orders covered by the l ibrary . 

We used the Yonsei-Yale (YY) isochronesS ()200ll ) to 
determine the physical properties of the host star using 
the stellar parameters derived with SPC. The isochrone 
analysis made used of the stellar effective temperature 
Teff = (5576 ± 50) K, the surface gravity \og{g) = (4.38 ± 
0.10) (cgs) and the metallicity [Fe/H] = -0.28 ± 0.08 
from the SPC analysis. The surface gravity acts as a lu- 
minosity indicator for the star. Whilst a/i?* derived from 
light curve fitting of the transits can also be used as a lu- 
minosity indicator ( jSeager fc IVlallen-Orn elas 2003), this 
parameter can often vary significantly between planet- 
only versus planet-with-moon models. Therefore, using 
\og{g) makes no prior assumption about which model is 
correct. Using the JVIonte Carlo method, we generated 



random realizations of these three parameters and found 
matching isochrones and their corresponding physical pa- 
rameters. The physical parameter realizations are then 
used to produce a kernel density estimator to calculate 
the mode along with the 68.3% confidence intervals. This 
process leads us to estimate 



0.774j:!] °36 Mq and 



i?* = 0.833 



+0.061 
-0.039 



i?0, which may be compared to the B12 



values of = 1.03 M© and = 0.95 i?£ 
5.6.3. Planet-only fits 

When queried from IVIAST, the KIC effective temper- 
ature and surface gr avity were rep orted as Teff = 5497 K 
and logg = 4.497 (|Brown et all 120 lit ). Both of these 
estimates are consistent with our SPC determination 
and we decided to use these values to estimate the 
quadratic limb darkening coefficients, ui = 0.4848 and 
(ui + U2) = 0.6990. The initial two models we regressed 
were Vp and Vp.lDi where the former uses the theoretical 
limb darkening coefficients as fixed values and the latter 
allows the two coefficients to be free parameters. We find 
that logZ(Vp,LD) -logZ(Vp) = +237.56 ±0.25 indicat- 
ing that our limb darkening coefficients were not optimal. 
We thus chose to set the coefficients to the maximum a- 
posteriori values of ui = 0.1389 and {ui + 1*2) = 0.8147. 

KOI-303.01 has a period of Pp = (60.92884 ± 
0.00013) days (as determined by model Vp,ld) and ex- 
hibits 10 full transits from Q1-Q9 from epoch -6 to -1-6, 
except epochs 0, -f2 and -t-3 which are absent. As is 
typical for all cases, logZ(Jp) > logZ(Vp) indicating 
that allowing for 10 independent baseline parameters is 
unnecessary relative to a single baseline term. 

We find no evidence for TTVs in KOI-303.01, with 
logZ(J"TTv) - logZ(J'p) = -24.20 ± 0.18, which is for- 
mally an 6.6 cr preference for a static model over a TTV 
model. The timing precision on the 10 transits ranged 
from 1.9 to 3.0 minutes (from model J-"ttv ) and yields 
a flat TTV proflle, as shown in Fig. 2(f) We calcu- 



late a standard deviation of (5ttv = 3.8 minutes and 
Xttv ~ -'^^•^ '^'-"^ degrees of freedom. 

The TTV+TDV model fit, Vy, finds consistent tran- 
sit times with those derived by model Jt'tv- The tim- 
ings show no cl ear pattern or excessive scatter, visi- 
ble in Fig. 2(f)[ We therefore conclude that there is 
presently no evidence for TTVs or TDVs for KOT303.01. 
The standard deviation of the TDVs is found to be 
(^TDV = 3.2 minutes and we determine Xxdv ~ ^^-i for 
10-1 degrees of freedom. 

5.6.4. Moon fits 

A planet-with-moon fit, J-s^ is slightly favoured rela- 
tive to a planet-only fit at A(logZ) = -1-12.34 ± 0.16, 
or 4.59 fj, meaning detection criterion Bl is met. The 
zero-mass moon model is slightly preferred though at 
AlogZ = 3.13 ±0.17. 

Despite failing B2, the fits yield broadly physical pa- 
rameters and the zero-mass moon preference could be in- 
dicative of a low signal-to-noise TTV/TDV data at this 
stage. Using = 0.774 M© and R„ = 0.833 i?© from 
our spectroscopic observations and SPC analysis, we de- 



termine Mp = 19.6+^5';4Me for Rp = 2.QUZo,m&R® 
and Ms = Q.Alt^^llM^ for Rs = 0.958l°:;5^^ i?®- The 
solution is also clearly not a close-binary with Ms/Mp = 
0.0301 

0.018 ^'^'^ three mutual transits driving the fit, as 
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TABLE 19 

Transit times and durations for KOI-303.01. The model used to calculate the supplied values is provided in parentheses next to each 

column heading. BJDuTC times offset by 2, 400, 000 days. 



Epoch T [BJDuTc] (-^TTv) 



TTV [mins] (J"ttv) t [BJDutc] (Vy) TTV [mins] (Vv) T 



s] (Vv) TDV [mins] (Vy) 



-6 


55006.3643loooi7 


-4.3 ±2.5 


55006. 3646+o;ooi9 


-3.6 ± 2.7 


-5 


55067.2957t|j;^E;i^ 


-0.6 ±2.2 


55067.2953to:ooi7 


-0.9 ±2.4 


-4 


55128.2232 q QQ^^g 


-2.6 ±2.5 


Ci'^l 0000 + 0.0019 


-2.6 ±2.8 


-3 


55189. 1563t°goi7 


+3.6 ±2.5 


55189.1562tH0i« 


+4.1 ± 2.3 


-2 


55250.0824+|j-gJj}^ 


-0.4 ± 2.5 


55250.0824ljj!j|]ig 


+0.3 + 2.5 


-1 


55311.0139«:rj 


+3.4 ± 2.7 


55311.0138in°i» 


+4.2 ± 2.6 


+1 


55432. 8720to;ooi8 


+4.1 ± 2.6 


55432.8719lo:oo2i 


+5.1 ± 3.0 


+4 


55615. 6569t°go" 


+1.7+1.9 


55615.6566to:ooi6 


+2.7 + 2.3 


+5 


55676. 5837^00021 


-1.2 ±3.0 


55676.5829lo:oo26 


-0.9 ±4.1 


+6 


55737.5080t|j;!]E!^? 


-7.8 ±3.0 


55737.5074to;oo23 


-7.0 ±3.5 



380.5l|;| 
379.2t^-J 
380.lt«J 
363.5t^;4 

373.lt57 
375.4t5-;8 



370.1 



+6.6 



365.6"* 



-6.1 
,+4.9 
-5.3 

375.6t^-^ 



366.6" 



8.8 



+4.6 + 3.3 
+3.9 + 2.8 
+4.3 + 3.3 
-4.0 + 2.7 
+0.9 + 2.9 
+2.0 + 2.9 
-0.6 + 3.2 
-2.9 + 2.6 
+2.1 + 4.5 
-2.4 + 4.1 



evident from Fig.|9l Finally, the mass ratio posterior very 
slightly peaks away from zero suggesting a low signal-to- 
noise hint of a TTV/TDV mass signal. In conclusion, we 
evaluate that detection criterion B3 is satisfied and B4 
is unclear. 

5.6.5. Predictive power of the moon model 

At this stage, we considered KOI-303.01 to be a poten- 
tial candidate and further data may confirm/reject the 
signal. We therefore detrended the PA data for KOI-303 
from QlO-13 covering 5 new transits, which had recently 
become available at the time of writing. Before we at- 
tempted to rc-fit the updated data set, we extrapolated 
the light curve model fit into QlO-13 for both the 
maximum a-posteriori solution and 50 randomly sam- 
pled solutions from the joint posteriors, which are shown 
in Fig. |9l The maximum a-posteriori moon solution can 
be compared to the detrended data, where we compute 
= 1295.21 for 836 data points. For comparison, we re- 
peated the process for the maximum a-posteriori planet- 
only model and compute = 1276.85. To provide some 
reference, the two best fits yield a oi 11688.64 and 
11746.60 over Ql-9 for J^s and respectively over 10344 
data points. 

This simple test suggests the moon model lacks pre- 
dictive power, which is a major concern for accepting 
the moon hypothesis. In order to investigate the uncer- 
tainty in the best-fit models, we extrapolated the transit 
light curve for 50 realizations where the system parame- 
ters are randomly drawn from the joint posteriors. The 
50 realizations from the model J^s are shown in Fig. [9l 
The distribution of the values for these fits yields 
= 1300 ± 21 whereas repeating the process for 50 pre- 
dictions from yields 1280 ± 25, which again suggests 
that the planetary model has superior predictive power. 
We therefore conclude that KOI-303.01 does not satisfy 
detection criterion F2. On this basis, we reject KOI- 
303.01 as a possible candidate and consider that there 
is no evidence to support the presence of an exomoon 
around this target with present data. 

The three mutual events fitted by the moon model 
mo rphologically resem ble starspot crossing events (e.g. 
see iRabus et al.l 120091 ) . We consider that this is the 
most likely source of false-positive for KOI-303.01 in 
light of the unmet detection criteria and the presence of 
flux variations consistent with rotational modulations. 
These modulations reveal a complex set of periodici- 
ties when processed with a Lomb-Scargle periodogram 



TABLE 20 

Bayesian evidences of various fits for KOI-303.01. A description 
of the different models can be found in ^4.1\ 



Model, M 



\ogZ{M) 



Ml - M2 
= log2(A4i) 



logZ{M2) 



Planet only 
Vp 

Vp,LD 
Vp.MAP 



fits... 

69125.46 + 0.18 
69363.02 + 0.18 
69365.13 + 0.18 
69439.45 + 0.11 



Yp,ld — Vp_ = 

Vp.MAP — Vp 
J7p _ Vp MAP 



(+237.56 + 0.25) 
= (+239.67 + 0.25) 
= (+74.32 + 0.21) 



Planet with timing variations fits. . . 

J^TTV 69415.25 ± 0.14 Tttv - = 

Vv 69218.80 ±0.27 Vv - Vp map 



(-24.20 ±0.18) 
= (-146.33 ±0.32) 



Planet with moon fits... 

69451.79 ±0.12 J's - J'p = (+12.34 ± 0.16) 

J^S.mo 69454.92 + 0.12 J-g.MO - -^s = (+3.13 ± 0.17) 

J-g^RO 69447.02 ± 0.16 Ts]ko - = ("4.77 ± 0.20) 



(see Fig. [T0)) . indicative of significant differential rota- 
tion. Analyzing Ql-13 (except Q12 which is incoherent 
with the other quarters), the highest power occurs at 
25.92 days and the second highest at 30.41 days. We es- 
timated spot anomaly timings of BJDutc 2455189.173, 
2455250.013 and 2455615.604. These timings are incon- 
sistent with a single spot rotating every 25.92 days but do 
give excellent agreement with a single spot rotating every 
30.41 days. We therefore consider that KOI-303 exhibits 
differential rotation and the transit crosses a spot which 
rotates every 30.41 days. 

One major effect of spots are such crossing events yield- 
ing a false Rs/Rp signal. To a lesser degree, spot cross- 
ings can perturb the transit profile and cause small false 
timing variations, giving a false Ms/Mp signal. How- 
ever, in general this latter distortion will be small and 
already our timing analysis suggests a lack of TTVs. In 
light of this, we use the J"s.Rn model for our final AIs/Mp 
posterior, shown in Fig. 4(f) The 95% quantile from this 
posterior constrains Ms/M p < 0.21 and the 3cr quantile 
constrains Ms/Mp < 0.33. 

5.6.6. Summary 

We find no compelling evidence for an exomoon around 
KOI-303.01 and estimate that Ms/Mp < 0.21 to 95% 
confidence. This assessment is based on the fact the sys- 
tem fails the basic detection criterion B2 and is marginal 
for B4. Further investigation reveals that the candidate 
also fails the follow-up criterion F2 (see §4.5p . 
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Fig. 9. — From left-to-right then top-to-bottom we show the chronological sequence of transits of KOI-303.01. The first 10 panels show 
the Ql-9 data and the maximum a-posteriori light curve fit of a planet-only model (gray line) and a moon model (red line). The last 5 
panels show the QlO-13 data with 50 extrapolations (50 different shadings used) of the moon model overlaid (parameters randomly drawn 
from the joint posteriors), which do not exhibit significant predictive power. 
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Fig. 10. — Lomb-Scargle periodogram of the out-of-transit nor- 
malized flux of KOI-303 from quarters 1-11 and quarter 13. We 
zoom in on 10-100 days where significant power is identified, ex- 
hibitting a complex set of periodicities indicative of differential ro- 
tation. The peak at 30.41 days is consistent with the timings of the 
transit anomalies. 



TABLE 21 

System parameters for KOI-303. 01 from model Vp,LD> except for 
Ms/Mp which is derived from model J^s.RO- 



Parameter 



Value 



Derived parameters 
Pp [days] 
TO [BJDUTC] 
Rp/R^, 



60.92884l»;~Sf3 



2455371.94026 



+0.00052 



U.UZOU*_g QQQ44 



b 

(a/R.) 
i [deg] 
pi, [gcm~^ 
f [hours] 

Ml 

(ui + U2) 



^+0.27 
-0.24 
:t + i.2 
-11.1 
1-1-0.20 
-0.31 
1 79+0-34 
^- ' =^-0.72 

-1+0.065 
-0.063 



0.363 
70.61 
89. 71^ 
1.79l| 
6.150"* 



Physical parameters 

[Rq] (SPC) 
R* [Rq] (SPC) 

Rp [R®] 
Ms/Mp 
Sttv [mins] 
(Stdv [mins] 



rs 774+0.036 
'*-0.017 
+0.061 
0.039 
+0.17 
0.13 



0.833: 
2.30j 

< 0.21 (95% confidence) 

< 3.6 (95% confidence) 

< 3.9 (95% confidence) 



5.7. KOI-1816 

5.7.1. Data selection 

After detrending with CoFiAM, the PA and PDC-IVIAP 
data were found to have a 1.1a and 1.9 tr confidence of 
autocorrelation on a 30 minute timescale respectively and 
therefore both were acceptable (< 3cr). In general, we 
always prefer using the raw data and so we opted for the 
PA data in all subsequent analysis of this system. We 
note that only long-cadence data was available for this 
system. 

5.7.2. Planet- only fits 

When queried from JVIAST, the KIC effective temper- 
ature and surface gravity were reporte d as T^s = 4230 K 
and \ogg = 4.387 (iBrown et al.ll20lil ). Using these val- 
ues, we estimated quadratic limb darkening coefficients 
iti = 0.6747 and (U1+M2) = 0.7511. The initial two mod- 



els we regressed were Vp and Vp^ld where the former uses 
the aforementioned limb darkening coefficients as fixed 
values and the latter allows the two coefficients to be free 
parameters. We find that logZ(VpxD) — log-2^(Vp) = 
—0.24 ±0.20 indicating that the theoretical limb darken- 
ing coefficients are adequate and will be adopted in all 
subsequent fits. 

KOI-1876.01 has a period of Pp = 82.53238 ± 
0.00071 days (as determined by model Vp^ld) and ex- 
hibits 8 transits from Q1-Q9. As is typical for all cases, 
logZ(J^p) > logZ(Vp) indicating that allowing for 8 in- 
dependent baseline parameters is unnecessary relative to 
a single baseline term. 

We find no evidence for TTVs in KOI-1876.01, with 
logZ(J-TTv) - logZ(J'p) = -17.27 ± 0.15, which is for- 
mally an 5.5(7 preference for a static model over a TTV 
model. The timing precision on the 8 transits r anged 



from 6.0 to 7.9 minutes. The TTVs, shown in Fig. 2(g) 



show no clear pattern and exhibit a standard deviation 
of (5ttv — 10.8 minutes and Xxtv ~ ^^-^ ^'^^ ^"^ degrees 
of freedom. 

The TTV+TDV model fit, Vy, finds consistent tran- 
sit times with those derived by model J^ttv- No clear 
patter n or e xcessive scatter is visible in the data, shown 
in Fig 



2(g) 



We therefore conclude there is no evidence 
for TTWot TDVs for KOI-1876.01. The standard de- 
viation of the TDVs is found to be (5tdv = 14.2 minutes 
and we determine Xtdv = 4.6 for 8-1 degrees of freedom. 

5.7.3. Moon fits 

A planet- with-moon with, J^g , is preferable to a planet- 
only fit at a formally moderate significance level of 3.4 cr, 
close to our A a threshold of criterion Bl. Detection cri- 
terion B2 is certainly satisfied with J^s favored over the 
zero-mass moon model at 6.0 a and over the zero-radius 
moon model at 13.6 a (see Table [22|) . 

We may combine the posteriors from with the stel- 
lar parameters derived by B12 (Af* = 0.51 and 
i?^ = 0.49 i?®) to obtain physical parameters for the 
planet-moon candidate system. The planet's parameters 
are consistent with a dense Super-Earth/mini- Neptune 
with Mp = 34.2t^jMe for Rp = 2.08lt^;^^^ i?®. 
In contrast, the moon appears to have an unphysi- 
cally low-density with Ms = 0.13lo;JgMe for Rp = 
0.9951 

0.091 ^® ■ conclude that this moon fit fails de- 
tection criterion B3. The Ms/Mp posterior, shown in 



Figure 4(g) also fails to converge away from zero, failing 
criterion b4. The 95% and 3 u quantiles of this posterior 
are Ms/Mp < 0.012 and Ms/Mp < 0.022. 

The zero-peaked moon mass posterior suggests that 
the spurious signal is driven by the spurious moon's ra- 
dius. Inspection of the maximum a-posteriori realization 
of model J-s (Fig. Ilip reveals auxiliary transits, which 
cannot be due to starspot crossings. For this reason, 
similar to KOI-1857.01, we decided to investigate the 
predictive power of the moon model, to ensure we are 
not overlooking a physical signal. 

5.7.4. Predictive power of the moon model 

We detrended the PA data for KOI-1876 from QlO-13 
covering 4 new transits, which had recently become avail- 
able at the time of writing. We extrapolated the light 
curve model fit J-s into QlO-13 for both the maximum 
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Fig. 11. — From left-to-right then top-to-bottom we show the chronological sequence of transits of KOI-1876.01. The first 8 panels show 
the Ql-9 data and the maximum a-posteriori light curve fit of a planet-only model (gray line) and a moon model (red line). The last 4 
panels show the QlO-13 data with 50 extrapolations (50 different shadings used) of the moon model overlaid (parameters randomly drawn 
from the joint posteriors), which do not exhibit significant predictive power. 



a-posteriori solution and 50 randomly sampled solutions 
from the joint posteriors, which are shown in Fig. 1111 The 
maximum a-posteriori moon solution can be compared 
to the detrended data, where we compute = 621.9 
for 516 data points. For comparison, we repeated the 
process for the maximum a-posteriori planet-only model 
and compute — 577.7 i.e. substantially worse. 

This comparison suggests that the moon model has no 
significant predictive power, which is a major concern 
for accepting the moon hypothesis. We should expect a 
moon model to give a better by around > 25 based 
upon the improvement of the best fits on the Ql-9 data. 
For Ql-9, we found Ax^ =52.9 between the two best fits 
for 1129 data points. The 50 realizations of light curve 
predictions, shown in Fig. 1111 do not seem to show a 



convincing agreement with the data cither. Indeed, the 
distribution of the values for these fits yields x^ = 
631 ± 25 whereas repeating the process for 50 predictions 
from J^p yields 593 ± 20, which again suggests that the 
planetary model actually has slightly better predictive 
power. 

The lack of predictive power of model J-g (detection cri- 
terion F2) confirms that the moon model is spurious and 
we conclude the most probable origin of the spurious sig- 
nal to be time-correlated noise, similar to KOT1857.01. 
There is therefore no evidence for an exomoon around 
KOI-1876.01 based on the presently available data. Our 
final system parameters for this system are provided in 
Table [231 
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TABLE 22 

Bayesian evidences of various fits for KOI- 1 876. 01 . A 
description of the different models can be found in ^4-l\ 



Model, M logZ(M) 



Ml - M2 

= logZiMi) - \ogZ{M2) 



Planet only 
Vp 

Vp,LD 



fits. . . 

6956.16 ± 0.14 
6955.92 ± 0.14 
7003.43 ± 0.09 



Vp.LD - Vp = (-0.24 ± 0.20) 
J-p - Vp = (+47.27 ± 0.17) 



Planet with timing variations fits. . . 

J"ttv 6986.16 ± 0.12 J"ttv - .^p = (-17.27 ± 0.15) 

Vy 6874.39 ±0.21 Vy - Vp = (-81.77 ± 0.25) 



Planet with 

■^S,MO 
•^S,RO 



moon fits. . . 
7010.71 ± 0.11 
6916.01 ± 0.11 
6990.39 ± 0.12 



JFg - J-p = (+7.27 ±0.14) 
-^S,MO -^S = (-94.70 ±0.15) 
-^S,Ro --^S = (-20.32 ±0.16) 



TABLE 23 

System parameters for KOI-1876.01 from model Vp ldi except 
for Ms/Mp which is derived from model J-g. * indicates that a 
parameter was fixed. 



Parameter 



Value 



Derived parameters 
Pp [days] 
TO [BJDuTC 
Rp/R* 
b 

(a/R,) 

i [deg] 89.78 

p* [gcm~^] 4.5 

f [hours] 4.72 

ui 0.61 

(UI+M2) 0.64 



09 C-qOqO+O. 00073 

2455372. 2615logoi7 
U.U2384_(, 00034 

118t" 



35 

+0.16 
0.33 
+ 1.7 
3.0 
+0.17 
0.18 
+0.47 
0.37 
+0.24 
0.27 



Physical parameters. 



M. [Rq] 
R* [Rq] 
Rp [R®] 

Ms/Mp 
5xTy [mins] 
<5tdv [mins] 



0.51* 
0.49* 

r, r,r,r, + 0.229 
Z.ZUZ_o ggg 

< 0.012 (95% confidence) 

< 5.4 (95% confidence) 

< 30.8 (95% confidence) 



5.7.5. Summary 

We find no compelling evidence for an exomoon around 
KOI-1876.01 and estimate that Ms/Mp < 0.012 to 95% 
confidence. This assessment is based on the fact the sys- 
tem fails the basic detection criteria B3 and B4 as well 
as the follow-up criterion F2 (see §4.5p . 

6. DISCUSSION & CONCLUSIONS 

We have presented the first systematic search for the 
moons of extrasolar planets. In this work, we have fo- 
cussed on a sub-sample of seven transiting planet can- 
didates which our automatic target selection (TSA) al- 
gorithm identified as being not only dynamically viable 
satellite hosts, but also exhibiting sufficiently high signal- 
to-noise photometry that Earth-sized moons should be 
detectable. We also point out that this sub-sample was 
selected with the additional filters that the planetary 
candidates have radii < 6 and that no other plan- 
etary candidates are known to exist in the system. The 
former filter is applied since Jovian-sized candidates have 
a higher false-positive rate (|Santerne ct al...2012.) and the 



latter since the presence of perturbing planets compli- 
cates our analysis. 

We find no compelling evidence for an extrasolar moon 
in any of the seven systems. This determination is aided 
by the introduction of exomoon detection criteria, which 
we have proposed in this paper (see §4.5p . Although 
some perverse configurations of satellite orbits could still 
hide large moons in the systems analyzed, we are able 
to marginalize over the parameter volume to place esti- 
mated upper limits on the satellite-to-planet mass ratio, 
Ms/Mp (see Fig.El). 

In the cases of KOI-722.01, KOT1472.01, KOI-1857.01 
and KOI-1876.01, we estimate tight constraints with a 
95% upper quantiles of Ms/Mp < 0.04. Since all four 
objects are relatively low-radii (2.0 i?©, 3.9 i?©, 2.2 i?© 
& 2.2 i?0 respectively), these limits likely probe down to 
sub-Earth masses, although an exact determination can- 
not be made without an estimate of the planetary masses. 
For KOI-365.01, KOI-174.01 and KOI-303.01 we are un- 
able to derive tight upper limits due to the fits favoring 
close-binary type solutions. A satellite in close proxim- 
ity shows negligible transit timing effects since the orbital 
period of the moon is comparable to that of the transit 
duration. For this reason, almost any Ms/Mp value is 
allowed and our upper limits on this term for these three 
candidates are less useful. In the case of KOT303.01, 
this spurious signal is probably due to the presence of 
three starspot crossing events driving a spurious fit. In 
the cases of KOI-365.01 and KOI-174.01, the origin is 
less clear but we suggest here the possibility of erroneous 
limb darkening parameters, time-correlated noise or stel- 
lar activity inducing slight transit profile distortions. 

In all seven cases, we are able to derive upper limits on 
the presence of moon-induced transit timing variations 
(TTV) and transit duration variations (TDV) , which are 
typically of the order of a few minutes. We find that 
all of the candidates favor a static model over a model 
including TTVs. Transit times and durations are also 
made available and are visible in Fig. [2j Revised transit 
parameters are presented for all planetary candidates, in- 
cluding revised stellar parameters for KOI-303.01 based 
upon new observations using the ARCES spectrograph 
on the 3.5 m Astrophysical Research Consortium Tele- 
scope at the Apache Point Observatory. The light curve 
derived stellar densities from our fits are consistent with 
the KIC stellar classifications, suggesting no obvious ev- 
idence for a blend or a wildly eccentric orbit in any case. 

With only seven candidates analyzed and four yielding 
strong upper limits, we emphasize that our HEK survey 
is only just beginning and it is too early to draw any 
meaningful statistical conclusions at this stage. The four 
candidates with strong limits have radii between 2 i?© 
and 4i?0 and may be classed as Super-Earths/lVIini- 
Neptunes. Our results therefore cautiously suggest the 
preliminary conclusion that such objects do not acquire 
large moons with a high frequency. This may be be- 
cause the objects never acquired a large moon or alterna- 
tively because such large moons arc lost quicker than the 
timescales estimated in Barnes & O'Brien (2002), which 
we used for our target selection procedure. The latter 
argument may be investigated by inspecting larger radii 
planetary candidates since the hosts will be presumably 
more massive and thus capable of maintaining a large 
satellite for a greater duration. The relative paucity of 



II. The Hunt for Exomoons 



33 



Jovian-sized objects in the Kepler-sample (B12) means 
that they are generally less favorable for a HEK analysis 
(e.g. fewer bright host stars), but nevertheless we will 
survey these candidates in future work. 

The other unique characteristic of this survey was that 
we focussed on planetary candidates in systems where 
no other planetary candidates had been identified. This 
choice simplifies our TTV analysis but also represents a 
sample bias. Closely-packed systems with multiple plan- 
ets may translate to a greater possibility of large moons, 
perhaps due to a larger initial reservoir of planetessi- 
mals from which t o form such bodies . In contrast, the 
system Kepler-36 (jCarter et al.l I2012D offers a specula- 
tive counter-example since the only two known transit- 
ing bodies reside in 6:7 resonance and may have obtained 
such a configuration through the stripping of one of them 
as a primordial moon. Since only two transiting bodies 
are known, these objects would have appeared as a sin- 
gle transiting binary in a previous epoch, similar to the 
candidates focused on in this survey. It is clear that we 
are only just beginning to unearth evidence for or against 
large moons and future HEK surveys will test these hy- 
potheses by sampling different systems. 
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